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Abstract: One-instanton contributions to the correlation functions of two gauge- invariant 
single-trace operators in J/ = 4 SU(A^) Yang-Mills theory are studied in semi-classical 
approximation in the BMN limit. The most straightforward examples involve operators 
with four bosonic impurities (whereas examples with two-impurity operators pose technical 
problems). The explicit form for the correlation functions, which determine the anomalous 
dimensions, follows after integration over the large number of bosonic and fermionic moduli. 
Our results demonstrate that the instanton contributions scale appropriately in the BMN 
limit. We find impressive agreement with the .D-instanton contributions to mass matrix 
elements of the dual plane-wave IIB superstring theory, obtained in a previous paper. Not 
only does the dependence on the scaled coupling constants match, but the dependence on 
the mode numbers of the states is also in striking agreement. 
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1. Introduction 

According to [1] there is a very interesting limit of the AdS/CFT correspondence that 
relates a special sector of the jV=A supersymmetric Yang-Mills (SYM) theory to type IIB 
string theory in a maximally supersymmetric plane-wave background. A notable feature 
of this proposal is that it provides the first example of a gauge/gravity duality which can 
be studied in a quantitative way beyond the supergravity approximation. This is possible 
because string theory in the relevant background, which is obtained as a Penrose limit of 
AdSs x S 5 [2], can be quantised in the light cone gauge [3,4]. Moreover there exists a 
regime in which both the string and the gauge theory are weakly coupled, although there 
are subtleties associated with the different way in which the limits leading to this regime 
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are taken on the two sides. This has allowed very precise comparisons between perturbative 
corrections in the two theories [5,6]. 

The duality relates the string mass spectrum to the spectrum of scaling dimensions of 
gauge theory operators in the so called BMN sector of =/K=4 SYM. This consists of gauge 
invariant operators of large conformal dimension, A, and large charge, J, with respect to 
a U(l) subgroup of the SU(4) R-symmetry group. The duality involves the double limit 
A — > oo, J — > oo. The combination A — J, which is kept finite, is related to the string 
theory hamiltonian, 

A-J=-H {2 \ (1.1) 

where p is the background value of the R&R five- form and is related to the mass parameter, 
to, which appears in the light cone string action by to = pp-ol [3,4], where p- is the light 
cone momentum. 

The correspondence between the spectra of the two theories is thus the statement 



that the eigenvalues of the operators on the two sides of the equality (1.1) coincide. A 
quantitative comparison is possible if one considers the large N limit in the gauge theory 
focusing on operators in the BMN sector. As a result of combining the large N limit with 
the limit of large A and J, new effective parameters arise [7,8], which are related to the 
ordinary 't Hooft parameters, A and 1/N, by a rescaling, 

a 2 N J 2 
A' = ^, 92 = ^- (1-2) 

The correspondence relates these effective gauge theory couplings to string theory param- 
eters in the plane-wave background, 

m 2 = (pp^a) 2 = — , 4irg s m 2 = g 2 . (1.3) 

The double scaling limit, N —* oo, J — ► oo with J 2 /N fixed, connects the weak coupling 
regime of the gauge theory to string theory at small g s and large to. 

In this limit the leading perturbative corrections to the scaling dimensions of BMN 
operators have been successfully compared to the leading quantum corrections to the masses 
of the dual plane- wave string states [5-11], see also the reviews [12] for further references. In 
the present paper we will study one-instanton effects in the jV=A Yang-Mills theory. These 
will be compared with D-instanton [13] induced corrections to the plane-wave string mass 
spectrum that were computed in [14] in order to check the validity of the BMN conjecture 
in non-perturbative sectors. In the original formulation of the AdS/CFT correspondence 
very good agreement was found between the effects of instantons in the JY=A Yang-Mills 
theory and of .D-instantons in type IIB string theory in AdSs x 5 5 [15-17]. It is therefore 
of interest to see if similar agreement can be established in the BMN limit and whether 
the results of [14] can be reproduced from the study of instanton contributions to the 
anomalous dimensions of BMN operators. 

The possibility of testing the correspondence at the non-perturbative level is especially 
relevant since several aspects of the perturbative tests of the duality are only partially 
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understood. A precise holographic formulation of the duality connecting the dynamics of 
the two theories beyond the identification of the spectra is still lacking and even the explicit 
tests at the level of the spectrum are not comprehensive. A limited class of states/operators 
has been studied and agreement has been explicitly verified only at leading order in g 2 (the 
planar limit). This is the limit in which the string is free and the A' perturbation series on 
the gauge side reproduces the free string spectrum. The first non-planar contributions, of 
order X'g^ or one-loop on the string side, have also been compared successfully (although in 
this case simplifying assumptions were made about the one-loop string theory calculation, 
which have since been questioned [18]). The systematics of the perturbative expansion 
beyond these leading order contributions has not been studied and the fact that the double 
scaling parameters (|L^) that arise at low orders indeed represent the correct expansion 
parameters at all orders remains a conjecture. 

Results obtained in different but related limits of the AdS/CFT duality, both in string 
theory [19] and on the gauge theory side [20], suggest the possibility that BMN scaling, 
i.e. the order by order reorganisation of the perturbative expansion into a double series in 
A' and g 2 , might break down at higher orders. In the strict BMN sector the scaling (|1.2j ) 
has been verified to three loops in perturbation theory [21], but a deviation was observed 
in a related matrix model calculation [22] at four loops. 

In this paper we will show that instanton contributions to the conformal dimensions 
of BMN operators display BMN scaling. Two-point functions of BMN operators computed 
in the semi-classical approximation will be shown to be in striking agreement with the 
D-instanton induced two-point amplitudes computed in [14]. The agreement includes not 
only the dependence on the parameters A' and g 2 , but also the dependence on the mode 
numbers characterising the states. In particular the agreement with [14] in the mode 
number dependence is highly non-trivial and requires dramatic cancellations. These results 
combined with the three loop perturbative result provide substantial evidence indicating 
that BMN scaling should persist at all orders. 

Instanton contributions to the anomalous dimensions of BMN operators are extracted 
from two-point correlation functions computed in the semi-classical approximation. Con- 
formal invariance determines the form of two-point functions of primary operators, & and 
&, to be 

^( Xl )ff(x 2 )) = _ C ^ )2A , (1.4) 

where A is the scaling dimension. In general in the quantum theory A acquires an anoma- 
lous term, A(g YM ) = Aq + 7(<7 Y m)- At weak coupling the anomalous dimension "f(g YM ) is 
small and substituting in (|1.4[) gives 



xxW{x 2 )) = — _ Ao (1 - 7 (<? YM ) log [A (xi - x 2 f] +■■■), (1.5) 

where A is an arbitrary renormalisation scale. As a function of the coupling constant the 
anomalous dimension admits an expansion consisting of a perturbative series plus non- 
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perturbative corrections. The generic two-point function at weak coupling takes the form 

<^l)^2)) = {X1 C{9Y X ^2A " 4m7 W [A 2 (X! - X 2 f] 

+ ...- e 2— 7 (inst) lQg [ A 2 (xi _ X2) 2] + . . . \ (L6) 



Therefore perturbative and instanton contributions to the anomalous dimension are ex- 
tracted from the coefficients of the logarithmically divergent terms in a two-point function. 
When there is more than one operator with the same quantum numbers operator mixing 
occurs. In this case the resulting set of two-point functions determines a matrix of anoma- 
lous dimensions and the eigenvalues of this matrix are the physical anomalous dimensions. 
The issue of operator mixing was first discussed in the context of the BM N limit in [23] . 

The procedure for calculating the instanton-induced contribution to the anomalous 
dimensions in semi-classical approximation is as follows. The gauge-invariant operators 
in the BMN sector are defined by colour traces involving a large number of elementary 
scalar fields together with a finite number of bosonic or fermionic 'impurities'. In the semi- 
classical approximation correlation functions of such operators are computed by replacing 
the fields by the solution to the corresponding field equations in the presence of an instan- 
ton, expressed in terms of the fermionic and bosonic moduli, and integrating the resulting 
profiles over these moduli. These moduli encode the broken super confer mal symmetries 
together with the (super) symmetries associated with the orientation of a SU(2) instanton 
within SU(iV). For large N integration over these moduli is carried out by a saddle point 
procedure (as in [24]). 

The general structure of the anomalous dimensions of gauge invariant operators in the 
t /K=4 Yang-Mills theory with SU(iV) gauge group is an expansion of the form 

oo oo 

7 ( 5ym ,e,N) = j2 < e * (ao oil + E E Kf } w a 2 z * 2mrK + c - c - ], a- 7 ) 

ri=l K>0m=0 

where r = 4- + i^ 2 -. The double series in the second term in (|1.7j) contains the contribu- 
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YM 

tions of multi-instanton sectors as well as the perturbative fluctuations in each such sector. 
One reason for studying instanton effects even though they are exponentially suppressed 
in the small coupling limit, is that they determine the dependence on the #-angle in jV=A 
SYM. They therefore play an essential role in implementing S'-duality which is a symmetry 
of the theory, just as D-instantons are crucial for the S'-duality in type IIB string theory. 

If the BMN sector of the gauge theory scales appropriately (|F7]) becomes a series in 
the scaled couplings A' and g 2 - In particular, we will show that the leading one-instanton 
contribution to the two-point functions of a class of four impurity BMN operators scales as 
it should in the BMN limit and has the form (l/ninj) 2 ^ exp (—8ir 2 /g2^' where 
the integers n\ and n 2 correspond to the mode numbers of the dual string state. This 
result is in striking agreement with the corresponding D-instanton induced mass matrix 
on the string side found in [14] . 

For certain other classes of operators the leading one-instanton contribution vanishes 
and the first non-zero correction is of higher order in A' (or a lower power of m in the 
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string calculation). In such cases the calculation requires knowledge of a non-leading term 
in the scalar solution - a term involving six fermionic moduli (whereas the leading term is 
quadratic in fermionic moduli) . We have not evaluated the precise form of this contribution 
and so have not determined the precise form of the matrix elements in these cases. However, 
there is strong evidence that these also match the string calculations. For example, for two 
impurity operators, with some mild assumptions about the manner in which the fermion 
moduli are distributed in the profile of the operators, we will find a contribution to the 
two-point function of the form A' 2 ^ 2 exp (— 8ir 2 /g2\' + i&) , in accord with expectations 
from the string side. Later we will comment on the systematics of the expansion in the 
one-instanton sector and on how the higher order corrections can give rise to a double series 
in A' and 52- 

This paper is organised as follows. In section || we review some general aspects of the 
jV=A: Yang-Mills theory and the BMN limit. The general method for evaluating instanton- 
induced contributions to two-point functions of BMN operators in terms of zero modes and 
the integration over super-moduli is described in section |3| The manner in which the profiles 
of the fields depend on these moduli is presented in section |3|. In section |5| we consider 
some specific examples of two-point functions of BMN operators and derive expressions 
for the anomalous dimensions that arise after integration over all the moduli. We first 
consider the case of two-impurity operators (which presents the technical difficulty alluded 
to above) and then four-impurity operators. We conclude with a discussion in section |(| 
which includes a comparison with the string results in [14]. Some technical details of the 
calculations are presented in the appendices. 

2. Fields and operators in the BMN limit 

The purpose of this section is mainly to present the notation used in the paper and to 
define the dictionary to be used for the comparison with string theory in the plane-wave 
background. We will only consider a small set of BMN operators with scalar impurities 
which are dual to the string states studied in [14]. A more detailed discussion of the 
various types of operators relevant for the comparison with string theory in the plane-wave 
background can be found in the review papers [12]. 

2.1 Fields in ^=4 SYM 

The ,jV=A multiplet comprises six real scalars, i = 1, ... ,6, four Weyl fermions, A^, 
A = 1,... ,4, and a vector, A^, with field strength F^, all transforming in the adjoint 
representation of the gauge group. These are the building blocks used to construct gauge 
invariant composite operators which are classified according to the irreducible representa- 
tions of the superconformal group, SU(2,2|4). The latter are identified by the quantum 
numbers (A,ji, j%\ a, b, c) of the maximal bosonic subgroup SO(2,4)xSO(6), where A is 
the scaling dimension, j\ and j'2 the Lorentz spins and [a, b, c] the SU(4)~SO(6) Dynkin 
labels. 

Under the SU(4) R-symmetry group the scalars transform in the 6, the fermions in the 
4 (and their conjugates in the 4) and the gauge field is a singlet. It is often convenient to 
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label the scalars by an antisymmetric pair of indices in the 4, ^ AB \ subject to the reality 
condition 

<PAB = {v AB Y = \e A BCDV CD - (2.1) 
The two parametrisations of the jV=^ scalars, (p % and (p , are related by 

,%i 1 ,„AB ,„AB 1 sr^AB-i in o\ 

= —^^AB'P > <P = -j= S i V > ( 2 - 2 ) 

where (S^" 8 ) are Clebsch-Gordan coefficients projecting the product of two 4's (4's) 
onto the 6. They are defined in appendix [A]. The representation of the scalars in terms 
of the (p AB fields is the most convenient for instanton calculations since, as we shall see, it 
makes manifest which fermion zero modes in a correlation function can be soaked up by 
each scalar field. 

In the limit relevant for the comparison with string theory in the plane-wave back- 
ground the symmetry group is a contraction of the original group and the operators are 
classified according to representations of the bosonic subgroup SO(4)xSO(4)xU(l)xU(l). 
We shall denote by V the dilation operator and by J the U(l) generator selected by the 
Penrose limit in the dual AdS background. The SO(4)xSO(4)xU(l)xU(l) quantum num- 
bers are (si, S2; s' 1; s' 2 ; A, J), where A and J refer to V and J and the spins Sj and s\ refer to 
the two SO (4) factors. These can be considered to be respectively subgroups of the original 
SO(6) and SO(2,4) groups. This identification is not completely correct. The generators 
of the two SO(4)'s corresponding to the isometries of the dual string background, Gi, are 
related to the generators of the Euclidean Lorentz group and to those of an SO (4) subgroup 
of the R-symmetry group, Gi, i = 1, 2, by a similarity transformation, Gi = TGiT^ 1 . This 
distinction, however, will not be relevant for our analysis. 

Since a precise formulation of the gauge theory dual to the plane wave string theory is 
not known, the rules for the decomposition of the jY=\ fields according to representations 
of SO(4)xSO(4)xU(l)xU(l) are determined by the quantum numbers of the dual string 
excitations. The gauge invariant operators corresponding to states in the string spectrum 
will be discussed in the next subsection. String excitations created by bosonic and fermionic 
oscillators are associated respectively with the insertion of bosonic and fermionic elementary 
fields ("impurities") in composite operators. 

Bosonic excitations in the plane wave string theory originate from the vector of SO(8) 
which decomposes under SO(4)xSO(4) as 



8* 



(2.3) 



i.e. they are vectors of one SO(4) and singlets of the second or vice versa. Correspondingly 
in the jV=^ theory the six real scalars are reorganised into a complex field, Z , and its 
conjugate, Z, which are singlets of SO(4)xSO(4) and have A = 1 and J = ±1 respectively, 
and four real fields which transform in the 4 = (5,5) of the first SO (4) and are singlets 
with respect to the second and have J = and A = 1. The insertion of the four real 
scalars in a composite operator corresponds to the insertion of bosonic creation operators 



- 6 - 



with an index in one of the two SO(4) factors in the dual string state. States created by 
bosonic oscillators which are vectors of the second SO(4) correspond to operators involving 
insertions of D^Z. The fields D^Z are in the 4 = (i, i) of the second SO(4) and have 
J = 1 and A = 2. Explicitly the scalar fields are 



z = 


i> x - 


= , 


Z = 


A- 


= 2<p 23 












1 


v x = 




V 2 = 
















i 


<. 3 = 






V 6 


= V2 



(<. 12 - <. 34 ) , 



(2.4) 



Here, for convenience of notation, we have introduced the scalars (p l , i = 1, . . . , 4, related to 
four of the (p 1 's by a relabelling. This should not cause any confusion since in the following 
we shall only work with ( |2. , ) and we shall not need the SO(6) fields (A. 5). 

Unlike the scalar fields the fermions transform non trivially with respect to both 
SO(4)'s. The four Weyl fermions of the jV=^ SYM theory, A^, transform in the 4 of 
SU(4)~SO(6) and their conjugates, A^, transform in the 4. Their decomposition with re- 
spect to SO(4)xSO(4)xU(l) is dictated by the SO(4)xSO(4) decomposition of the SO(8) 
fermions of the light-cone string. The type IIB fermions transform in the 8 S of SO (8), 
which under SO(4)xSO(4) decomposes as 



8, 



1 

°'2 



1 

2'° 



5'° 



(2.5) 



In terms of the 8 S fermions S a and S a the decomposition is achieved via a projector [3], 

i(i±n), 

s a ^(s-y a e(s + )i. (2.6) 



The Yang-Mills fermions, A^f, have A 
A 



3 and U(l) charge A for A 



1, 4 and 



,,•"■'"*—» <«i.m « , . , v...,,, S v 2 . . ..-..< him 2 for 

= 2,3. Similarly their conjugates, A^, have A = § and U(l) charge \ for A = 2,3 and 
— | for A = 1,4. To match the quantum numbers of the string oscillators we choose the 
following decomposition 



K^^'J ©V> i 



1,4, 



(2.7) 



where the fermions ip^a are defined as 



(M + A)_i 



(2.8) 



where the matrix M + is related to the matrix II used in the plane-wave string theory to 
project the SO(8) fermions onto chiral SO(4) xSO(4) spinors. The spinors V' - ? and \ 

transform under SO(4)xSO(4) in the (2 L ;2 L ) = [(±,0) ; (±,0)] and have (J 
and (J 



|) respectively. 
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Similarly 

Ai^t?.,eCr, « = 2 >3 (2.9) 

where 

4,-J h = (M-\y\_ (2.10) 

2 2 

and M~ is also related to the projector used to define the SO(4)xSO(4) spinors in the 
dual plane-wave string theory. The fermions . and ^)~\ a transform in the (2^; 2^) = 

+ 2' a — 2 

[(0, ^) ; (0, ^)] representation and have respectively (J = \, A = §) and (J = — |, A = |). 
Some aspects of instanton contributions to BMN operators involving fermionic impurities 
will be discussed in [25]. 

The field strength, F^ u , is a singlet with respect to the first SO(4) and decomposes 
into F^ v transforming in the 3~ = (1,0) and 3 + = (0, 1) with respect to the second. F^ v 
both have J = and A = 2. 

In the string amplitudes in the plane- wave background P+ and P_ are conserved, so the 
operators in the gauge theory are conveniently classified according to the dual quantities, 
i.e. A — J and A + J respectively (A + J is actually infinite in the limit; it is proportional 
to but the proportionality constant diverges). Table [l] summarises the A, J and 
SO(4)xSO(4) quantum numbers for the ,yV=4 elementary fields. The notation SO(4)r 
and SO(4)<7 has been introduced to denote the SO(4) groups descending from the original 
SO(6) (R-symmetry) and SO(2,4) (conformal) groups respectively. 



Field 


A 


J 


A - J 


A + J 


SO(4) R 


SO(4) c 


Z 


1 


1 





2 


(0,0) 


(0,0) 


Z 


1 


-1 


2 





(0,0) 


(0,0) 




1 





1 


1 


(2' 2) 


(0,0) 


D^Z 


2 


1 


1 


3 


(0,0) 


(2' 2) 




3 
2 


l 

2 


1 


2 


(io) 


(i.o) 




3 
2 


1 

2 


2 


1 








3 
2 


1 

2 


1 


2 


(0,1) 


(0,1) 




3 
2 


1 

2 


2 


1 


(°,*) 


(0,2-) 




2 





2 


2 


(0,0) 


(1,0) 


F + 


2 





2 


2 


(0,0) 


(0,1) 



Table 1: SO(4)xSO(4)xU(l)xU(l) quantum numbers of the ,jV =4 elementary fields. 
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2.2 BMN operators 

The composite operators dual to states in the spectrum of string theory in the plane wave 
background are also classified in terms of the same quantum numbers. In particular, A — J, 
which is dual to the light-cone hamiltonian, measures the number of "impurities" and will 
be used to classify the operators 

At finite J and A the selection rules of the Jf=k theory, implied by the superconformal 
symmetry, apply. So only two-point functions of (primary) operators of the same dimension 
can be non-zero. More precisely the two-point functions that are relevant for the calculation 
of anomalous dimensions are 

<^,A>)< )A >)> , (2-11) 

where the subscripts denote the SU(4) representation and the scaling dimension. The 
SU(2,2|4) symmetry imposes A,; = Aj and = rj so that both A and J are conserved in 
two-point functions. In the case of the U(l) charge J this means that the two operators in 
a non-zero two-point function must have equal and opposite charges. 

Gauge invariant composite operators which correspond to physical string states are of 
the form 

1 J 

ffh—ik _ 2m(p 1 n 1 +p2n 2 A \-p k n k )/J 

J;rn...n k i — r 



Pl<P2< — <Pk 

Tr (Z Pl X ei Z P2 ~ Pl X e2 ■ ■ ■ z Pk ~ Pk -^X lk Z J - pk ) 
J 

2-7ri[(n2H \-n k )q2+(nz-\ hnfc)iJ3H \-n k q k ]/J 



/ „2 7\T\ J+k 

J A- .7-1 ( YM \ q 2 ,-,q k =0 



C 



/ <?2H — \-qk<J 
xTr (z J - {q2+ - +qk) X ei Z q2 Xe 2 - ■ ■ Z Qk X^ , (2.12) 

where A — J denotes the total number of impurities. Here the cyclicity of the trace has 
been used and, after the change of variables, pi — > q%, pi — ► qi — qi-i (i = 2, . . . , k), the 
sum over q\ has been performed resulting in the condition 

n\ = -(n 2 H hrifc). (2-13) 

In ( |2.12| ) the X^s denote generic impurities, i.e. any of the elementary fields discussed in the 
previous subsection. String states dual to these operators are created acting on the vacuum 
with bosonic and fermionic oscillators. The integers ri\, . . . , n& in ( p. 12 ) are identified with 



the mode numbers in the dual string state and the relation ( 2.1 3| ) corresponds to the level 
matching condition obeyed by the physical string states. 

String creation operators are in one to one correspondence with A — J = 1 impurities, 
see table |]. Bosonic oscillators d l _ n and a^ n in the (4; 1) and (1;4) of SO(4)xSO(4) 
correspond to the insertion of ip l and D^Z impurities respectively . Fermionic oscillators, 



Here we are using a different notation with respect to [14], where the oscillators at n were denoted by 



-9- 



SZ n and S_ n , in the (2^; 2i) and (2/j; 2/j) correspond to the insertion of i[) f and ^jTi 



a 



2' 1 2 



impurities respectively. In string theory for each type of excitation one must consider left- 
and right-moving oscillators. These correspond to the insertion of the same field, but with 
the associated rn in the phase factor in ( 2.12| ) being respectively positive or negative. 



The normalisation of operators of the form ( p. 12 ) involving only A — J = 1 impurities 



is such that their tree level two-point functions are of order 1 in the BMN limit, N — > oo, 
J — > oo with J 2 /N finite. Operators involving A — J = 2 impurities have vanishing two- 
point functions at tree level because for equal total A — J they are normalised by the 
same prefactor but their definition involves fewer sums. Therefore these operators do not 
correspond to independent degrees of freedom in the BMN limit. In some cases, however, 
the insertion of A — J = 2 impurities is necessary in order to construct combinations which 
are well behaved in the double limit N — ► oo, J — > oo at higher orders in perturbation 
theory. For instance it is necessary to consider terms in which pairs of tp 1 impurities are 
replaced by a Z insertion in order to cancel divergences which arise at the level of the 
leading non planar perturbative corrections [9]. 

Operators with A — J = 0, 1 are protected and so their two-point functions do not 
receive instanton contributions. At the level of two and more impurities the situation is 
more interesting. The spectrum is significantly richer and more importantly there appear 
unprotected operators. In the following we shall only discuss a small selection of gauge 
invariant composite operators with scalar impurities which are dual to the string states 
analysed in [14]. A complete analysis would require computing the two-point functions 
involving all the operators in each sector and then diagonalising the resulting matrix of 
anomalous dimensions. We shall not carry out this program in this paper, but we shall 
concentrate on a few specific cases which illustrate the striking agreement with the correc- 
tions to the string mass spectrum calculated in [14]. The generic operator with k scalar 
impurities is of the form 

1 J 

fiii...i k _ 2iri[(mA \-n k _ 1 )p l +(n 2 -\ \-n k _ 1 )p 2 A \-n k _ 1 p k _ 1 }/ J 

Jk-l I "ym" \ Pi,-,Pfe-i=0 

^ 87r ' Pl+-+Pk-I<j 

xTr (z J - ( - Pl+ - +Pk - 1 ^ ll Z pl ip i2 ■ ■ ■ , (2.14) 

In our discussion we shall only consider operators with an even number of impurities. This is 
because operators with odd A — J, which receive perturbative corrections [26], are expected 
not to receive contributions in the one-instanton sector. This is a prediction following from 
the calculation of string amplitudes in [14]. In the plane- wave string theory the absence 
of instanton contributions to two-point amplitudes of states with an odd number of non- 
zero mode excitations is a straightforward consequence of the structure of the D-instanton 
boundary state. In the ^V=4 theory, however, the corresponding statement is far from 
obvious. 

2.2.1 Two impurity operators 

No field in the JV=4 multiplet has negative A — J hence the two impurity operators 
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are obtained with the insertion of either two A — J = 1 fields or of a single field with 
A — J = 2. Because of the normalisation (2.12) only operators with two A — J = 1 
insertions are relevant in the BMN limit. 

Even restricting the attention to SO(4)c singlets, already at the two impurity level 
there is a rather rich spectrum of operators, which becomes even richer when multi-trace 
operators with the same quantum numbers are taken into account. In the SO(4)r singlet 
sector one can construct gauge invariant operators in the representations (0, 0) = 1, (1, 0) = 
3 + , (0, 1) = 3 and (1, 1) = 9 of SO(4)#. The singlet can be realised with the insertion 
of two scalars, two gauge fields (through covariant derivatives) or two fermions of the 
same chirality. Operators in the 3^ include combinations of two scalar or two fermionic 
impurities. The 9 can only be obtained with the insertion of two scalar impurities. 

The operators with two tp 1 insertions are 2 



7 l;J;n 



7 3±;J;n 



9;J;n 



2 N \ J+2 
7 / YM 

J 1 "8^" 



£ e 27Tipn/J Tr (Z J -V^V ) " Tr (Z J+1 Z) 

p=0 



1 



J 



J2 e 2nipn/J r m ^ f Z J- Plp [k ZPlp l] 



J 



O AfX J+2 

YM \ P=0 

8tt 2 



(2.15) 



(2.16) 



(2.17) 



1 



2nipn/J 



Tr ( Z J - p v {i Z p ^ ) 



Tr (V-V^V) 



where the projectors onto the 3 + and 3~ are defined as r^ kl = \ (<5 lfc <5 J ' — S d 5^ k ± e 1 ^ 1 ). 
The singlet operator ( p.!5| ) provides an example of what mentioned earlier about Z in- 
sertions. In order to define a well behaved BMN operator it is necessary to consider the 
combination in (2.15). The second term is needed to cancel a divergent contribution to the 
two-point function of &\ j n arising at the leading non planar level [9]. 

The operators (gig)- (gig) are dual to string states in the plane-wave background of 
the form 



a % _ n a l _ n \Q) h 

rT l * k - n a l _ n \ 

a-n°?-n\ Q )h, 



(2.18) 
(2.19) 
(2.20) 



where |0)^ denotes the BMN ground state and the indices . . . are taken to be in one 
of the two SO(4) factors (to be identified with SO(4) R ). The integer n in (gl|)-(gl^) 

2 Here and in the following we use square brackets to denote antisymmetrisation, curly brackets to denote 
symmetrisation and subtraction of the trace and parentheses to indicate symmetrisation without subtraction 
of the trace part. 
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corresponds to the level of the dual string excitation. The n = operators are protected 
and correspond to supergravity states. 

As already remarked, in order to compute the instanton induced anomalous dimen- 
sions of the various operators in each sector it is in principle necessary to diagonalise the 
appropriate matrix. In the following we shall not consider the problem of mixing between 
single- and multi-trace operators, since it is a subleading effect in the large N limit. In 
general, however, at the instanton level mixing occurs at leading order among all the single 
trace operators in each sector [27]. In the case of the two impurity operators it has been 
shown [28] that all the operators in different sectors have the same anomalous dimension 
as a consequence of superconformal invariance. Therefore in the following we shall only 
analyse the single operator in the 9 for which there is no mixing to resolve. Superconformal 
symmetry guarantees that the results apply to operators in others sectors as well. 

In general the problem of resolving the operator mixing and computing anomalous 
dimensions can be drastically simplified using the constraints imposed by superconformal 
invariance, in particular, the fact that all the operators in a multiplet have the same 
anomalous dimension as the superconformal primary operator. 

2.2.2 Four impurity operators 

The number of independent BMN operators grows very rapidly with the number of im- 
purities and at the four impurity level the spectrum is already extremely rich. Bosonic 
operators in the SO(4)c* singlet sector exist in the following representations of SO(4)^ 

(0,0) = 1, (0,1) =3+, (1,0) =3-, (0,2) = 5+, (2,0) =5-, 

(1,1) = 9, (1,2) = 15+, (2,1) = 15", (2, 2) = 25. (2.21) 



Operators relevant in the BMN limit involve four A — J = 1 impurities. The combinations 
which contribute to SO(4)c scalars are listed (up to permutations of the four fields) in 
table | 3 . 



(i) ip 1 f 3 \p k if 1 


(**) t^^D^ZD^ZD^ZD^Z 


(Hi) ^ipW^ZD^Z 


(iv) ip l tpi\p- aa ip- b 


(v) ^ 3 ^t* 


(vi) D ll ZD»Zii~ aa ip- b 


(vii) D^ZDoZj+tft* 


(vizi) ^D^-^t" 


(ix) f ""fc'f 









Table 2: A — J = 4 combinations of impurities 



The SO(4)c singlet sector contains the largest number of operators, involving all 
the combinations (i)-(xi) in table Operators in the 3 + can contain (i) , (ii) , (iv)- 
(vi), (viii) and (ix) and those in the 3~ (i), (ii), (iv), (vii), (viii) and (x). The 9 involves 
(i), (ii), (iv), (v) and (viii). Operators in the 5 + and 5~ can be obtained from (i), (iv) and 

3 ^imm3M4 is a projector onto the singlet, i.e. <5 M1M2 5 M3M 4, <5 M1M3 <5 M2M4 or £ M1M2M3M4 . 
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(ix) and from (i), (v) and (x) respectively, those in the 15 + and 15 from (i) and (iv) and 
from (i) and (v) respectively. In the 25 there is only one operator corresponding to the 
combination (i) with indices fully symmetrised. 

In the following we shall concentrate on a few specific two-point functions correspond- 
ing to the amplitudes computed in [14]. This will be sufficient to show how instanton 
contributions to gauge theory correlation functions precisely reproduce various features 
observed in string theory amplitudes. The operators we study in detail are those contain- 
ing four scalar impurities. These are of the form 

t r - J 

ff _ \ r „27ri[(rai+n2+n3)pi + (ri2+n3)p2+"3P3]/ J 

" / r;J;ni,n 2 ,n 3 — y — — 7 e 

/ „2 fj\ J+4 — ' 

5 / YM \ Pi,P2,P3=0 

x ' pi+P2+P3<J 

xTr ^ J -( pi+P2+P3 V^ P V j ^ p V fe ^ P3 ^) > (2-22) 

where t^ kl is a projector onto the representation r of SO(4)#. In particular in the singlet 
sector there are three independent operators in this class, corresponding to the three tensors 

tijkl = e v kl ' l m = S V S kl . tfjkl = SikSjl ■ (2-23) 



In section 5.2 we study instanton contributions to two-point functions of operators of the 



type (2.22). We discuss in detail the case of the singlet corresponding the choice of the 



ijkl 

projector. We will show that the dependence on both the coupling constants, A' and g2, 
and the mode numbers, rtj, is in exact agreement with the results of [14]. We also find that 
for all the operators with four scalar impurities in sectors other than the singlet instanton 
contributions to the matrix of anomalous dimensions are suppressed by powers of A'. This 
result is also in agreement with the string prediction. 



3. Instanton contributions to two-point functions 

In this section we recall some general aspects of the calculation of instanton contributions 
to correlation functions, in particular two-point functions, in jV=A SYM. 

In the semi-classical limit correlation functions are evaluated using a saddle point 
approximation around the classical instanton configuration. In this limit the computation 
of expectation values reduces to an integration over the finite dimensional instanton moduli 
space as parametrised in the ADHM construction [29] . Before focusing on operators of large 
dimension, A, and R-charge, J, in the following sections, we briefly discuss the general 
formalism for extracting instanton contributions to the anomalous dimensions of gauge 
invariant local operators. Comprehensive reviews of instanton calculus in supersymmetric 
gauge theories can be found in [24,30] and instanton contributions to anomalous dimensions 
of scalar operators in j¥=^ SYM were studied in detail in [27]. 

Contributions to the (matrix of) anomalous dimensions are extracted from the loga- 
rithmically divergent terms in two-point functions. In the semi-classical approximation in 
the background of an instanton the two-point function of a generic local operator, 0{x), 
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and its conjugate takes the form 



(<@(x 1 )<ff(x 2 ))mst = J dfl- mst (m h ,tn { )e Sin8t G{x\\ m b) fflf )S'(x 2 ) m h , m f ) , (3.1) 

where we have denoted the bosonic and fermionic collective coordinates by and mf 



respectively. In (3J) d//i nst (mt,, mf ) is the integration measure on the instanton moduli 
space, Si ns t is the classical action evaluated on the instanton solution and G andG denote 
the classical expressions for the operators G and & computed in the instanton background. 

A one-instanton configuration in SU(iV) Yang-Mills theory is characterised by 4N 
bosonic collective coordinates. With a particular choice of parametrisation these bosonic 
moduli can be identified with the size, p, and position, xq, of the instanton as well as 
its global gauge orientation. The latter can be described by three angles identifying the 
iso-orientation of a SU(2) instanton and 4iV additional constrained variables, w U qi and w au 
(where u = 1, . . . , N is a colour index), in the coset SU(iV)/(SU(A r — 2)xU(l)) describing 
the embedding of the SU(2) configuration into S\J(N). In the one-instanton sector in the 
^V=A theory there are additionally 8N fermionic collective coordinates corresponding to 
zero modes of the Dirac operator in the background of an instanton. They comprise the 
16 moduli associated with Poincare and special supersymmetries broken by the instanton 
and denoted respectively by r\ A and £ aA (where A is an index in the fundamental of 
the SU(4) R-symmetry group) and 8iV additional parameters, v A and u Au , which can be 
considered as the fermionic superpartners of the gauge orientation parameters. The sixteen 
superconformal moduli are exact, i.e. they enter the expectation values (3.1) only through 



the classical profiles of the operators. The other fermion modes, v A and v Au , appear 
explicitly in the integration measure via the classical action, 5; ns t. This distinction plays a 
crucial role in the calculation of correlation functions. The v A and v Au modes satisfy the 
fermionic ADHM constraints 



w au v£ = 0, v Au w u6l = Q, (3.2) 

which effectively reduce their number to 8(N — 2). 

In the one-instanton sector the gauge-invariant measure on the instanton moduli space 
takes the form 

/ d/xphys e" 5inst (3.3) 

-AN n 4N „2vrir , 4 



where the instanton action is 



7T 2 



S inst = -2vrir + S iF = -2mr + -^^e ABC D& AB & CD (3.4) 

Z 9 YM P 

with 

gl M 2vr 2^2 
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In ( |3.3|) we have omitted an overall (AMndependent) numerical constant that will be rein- 
stated in the final expression. 

The two-point function (|3.1|) thus becomes 



{0(x x )0(x 2 )) 



-4N AN 2nir 

(N - l)l(N - 2)1 



( dpd 4 x fl dVd^d^-Vd 

A=l 



N-2-A p4N-13 



16q z p- 

e ym 



(xi;x ,p, rj, £, v, v) &(x 2 ]x , p, V, £, v, v) . (3.6) 



Following [16] the integration over the non-exact fermion modes can be reduced to a gaus- 
sian form introducing auxiliary bosonic coordinates, x\ * = 1 } • • • > 6, to rewrite the gauge 
invariant measure as 

i 



-4N 4N 2-nir 
11 ^YM 

(N — l)!(iV — 2)! 



dpd 4 x d 6 X \\d 2 V A d 2 i A d N ~ 2 v A d N - 2 y A 



4N-7 

p exp 



A=l 

-2p\Y + —XAB^ AB 



(3.7) 



where xab = ^^abX 1 an d the symbols T, l AB are defined in (A.l). 

The fermion modes v A and U Bu enter explicitly in the classical profiles of the operators 
in the instanton background as well as in the measure through the instanton action. It is 
thus convenient to construct a generating function as in [17], which allows to deal easily 
with the otherwise complicated combinatorics associated with the integration over v A and 
v Au . We introduce sources, d\ and $Au, coupled to v A and v Au and define 



Z [■&,■&]-■ 



4N 4N 2-Kir 
n ^YM 

(N - l)!(iV - 2)! 



dp d 4 X d 6 X 11 d V d 2 C A d N ~ 2 U A d"" V 



A=\ 



4N-7 

p exp 



9ym 



,B 



An 



Performing the gaussian integrals over v and v and introducing polar coordinates, 

6 



(r,n), 



2 2 

— r . 



i=l 



we find 



_ 2- 29 7r- 13 4 M e 2 ™ r 
Z ^= (Ar_i) !( ^ M -2)! 



d P d 4 x d 5 ^n d Vd 2 ^p 4iv - 7 



A=l 



(3-8) 



(3.9) 



(3.10) 



where all the numerical coefficients have been reinstated. In ( 3.10| ) we have introduced the 
density 



3r(0,0;n,r) =exp 



irr 



(3.11) 
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where the symplectic form VL AB is given by 

6 

n AB = tf B n\ J2( Qi ) 2 = 1 - ( 3 - 12 ) 

Gauge invariant operators depend on the v A and v Au variables via colour singlet bilinears. 
These arise in symmetric or anti-symmetric combinations transforming respectively in the 
10 and 6 dimensional representations of the SU(4) R-symmetry 

(sV^o ^ u< A u B ) = (z^Vf + u Bu i4) , (3.13) 
(u A u B ) 6 ^ 9^ = (u A -u B - v B "u A ) • (3-14) 

Using the generating function defined in ( |3.10| ) the v A v B bilinears in the operators & and 
G in ( |3.6| ) can be rewritten in terms of derivatives of t?; ft, r) with respect to the 
sources, $a and &b- The result for a two-point function in which the operator insertions 
contain a total of p {vv)q and q {pv)\a bilinears is of the form 

( ^ (2;i) ^ X2)) = (N-™KN-2)\ J d ^ d ^ II dVd 2 e A p 4iV - 7 (3.15) 



where & and & contain the dependence on the exact moduli, r\ A and £ , and on the bosonic 
collective coordinates after extracting the v A v B bilinears. Computing the r integral gives 

r 4 

(0(x 1 )0(x 2 ))~a( P ,q;N)g 8 Y +P + '>e 2mT / dpd 4 ^ d 5 U J] d V d 2 ^ / + ^ 5 

^ A=l 

f(U) ^(xnp^xoiri,^) 0(x2\p,xq;ti,£) , (3.16) 

where /(O) contains the dependence on the Q AB variables obtained from the derivatives 
of ft,r). The coefficient a(p,q;N) contains the N dependence and in the large N 

limit we find 

a{p,q;N) = — — ^ — tlL ( ]V p+ 2 + 0(N p+ 2 1 



(N - l)\(N - 2)1 

—ya+oa/AT),. (3.i7) 



From ( |3.16 ) and ( |3.17 ) it follows that the insertion of a (z^)io or (z/z^e bilinear in a 



correlation function produces a factor of <? YM or g YM y/N respectively. 

In computing the moduli space integrations in expressions for two-point functions of 



the type ( 3.16 ) it will prove convenient to calculate first the fermionic integrals over r\ A and 
£ A and the angular integration over the five-sphere. These give rise to selection rules that 
determine which operators receive instanton contributions to their scaling dimensions. In 



- 16 - 



particular since the superconformal modes are exact a correlation function can only receive 
instanton contribution if the operator expressions contain exactly sixteen such modes in 
the combination 

n m aA ) ■ (3.i8) 

A=l 

The integration over the five-sphere parametrised by the angular variables tt AB factorises 
and gives rise to further selection rules. It gives a non-vanishing result only if the SU(4) 
indices carried by the fi's in the two operators can be combined to form a SU(4) singlet. 
The SU(4) indices are originally carried by the fermion modes which are all in the 4, so the 
only possible singlet combinations correspond to products of s ABCD tensors. The generic 
five-sphere integral is of the form 

J d 5 n n AlBl . . . n A ^ B *" = c (n) ( £ AlBlA2B2 . . . e M ^ B ^ A ^ B ^ + permutations) , 

(3.19) 

where the normalisation constant c(n) is 

c (n = -i M. . (3.20) 

Equations ( |3. 18| ) and ( |3. 19| ) imply that a two-point function can receive a non-zero con- 
tribution only if the combined profiles of the two operators contain fermion modes of the 
four flavours with the same multiplicity. 

The bosonic integrations over the position and size of the instanton are left as a last 
step. In the case of two-point functions these integrals are logarithmically divergent, sig- 
nalling a contribution to the matrix of anomalous dimensions. 



4. Fermion zero modes 



In order to evaluate instanton induced correlation functions we need to integrate the clas- 
sical profiles of the relevant composite operators over the instanton moduli space. We 
are interested in the dependence on the collective coordinates and of particular relevance 
will be the way the fermionic modes enter into the expressions for the various fields. The 
zero-mode dependence in the elementary fields of the .JV=\ SYM multiplet was reviewed 
in detail in [27]. Here we briefly summarise the features which will be relevant for the 
analysis of two-point functions of BMN operators. 

The field equations of the t yf / =4 SYM theory admit a solution in which the gauge 
potential corresponds to a standard instanton of SU(iV) pure Yang-Mills theory and all 
the other fields vanish, 

A„ = 4, <p AB = \A = \*=Q. (4.1) 

However, the Dirac operator has zero modes in the background of this solution, i.e. the 
equation ]/)^ a X aA = has non-trivial solutions when the covariant derivative is evaluated in 
the background of an instanton. The general solution to the Dirac equation is linear in the 
instanton fermion zero modes. This non-trivial solution gives rise to a non-zero solution for 
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the scalar fields when plugged into the corresponding equation, D 2 ip AB = \^2[X A , X B ]. The 
latter admits a solution for the scalar which is bilinear in the fermion modes. Proceeding 
with this iterative solution of the field equations one generates a complete supermultiplet 
and further iterations give rise to additional terms with more fermion modes in each field. 
The general solution obtained through this procedure is schematically of the form 



n=0 n=0 
4n<8N 4n+2<8N 

^ = E A i 1+4nM > E ^ +4n) , (4-2) 

n=0 n=0 
4ra+l<8AT 4n+3<8iV 

where the notation $( ra ) is used to denote a term in the solution for the field <3? containing 
n fermion zero modes. It is also understood that in ( [4.2[ ) the number of superconformal 
modes in each field does not exceed 16 and the remaining modes are of v A and v type. 

In computing the expressions for gauge invariant composite operators we shall make 
use of the ADHM description in which the elementary fields are written as [N + 2] x [N + 2] 
matrices. In particular, the leading order term in the solution for the scalars ip AB is given 
explicitly in appendix The solution of the iterative equations becomes very involved 
after a few steps. However the flavour structure of the combination of fermion zero modes 
in each term can be determined without actually solving the equations and is sufficient to 
identify which operators can get an instanton correction to their scaling dimension. 

All the fermion zero modes, both the superconformal ones, r\ A and £ aA , and the modes 
of type v A and i> Au , transform in the 4 of SU(4). We shall denote a generic fermion mode by 
mf. The starting point for the construction of the instanton supermultiplet is the classical 
instanton, A^ , which has no fermions. The first term in A^ is linear in the fermion modes 

X$ A ~mf. (4.3) 

For the term 

^)AB in 

the scalar solution one finds 

,A2)AB _ [A B] fA .v 

(f y ~ m l m i , (4.4) 

i.e. the two fermion modes are antisymmetrised in order to form a combination in the 6. 
The 4 spinor contains three fermion modes in the combination 

\(3)<i BCD (a r\ 

A A ~ &ABCD m f m f m f > (4.t>J 

so that the component A^ has one mode of each flavour apart from A. Proceeding in the 
multiplet we find the quartic term in the solution for the vector, ^4^, which contains one 
fermion mode of each flavour in a singlet combination 

Aff ~ SABCD mfmfmfmf . (4.6) 

(5)A 

The following term is Aq , which has flavour structure 

~ £BCDE m { m f m { m f m f , (4.7) 
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i.e. it involves a mode of flavour A plus one of each flavour. Then we find <p^ AB that 
contains an antisymmetric combination of a mode of flavour A and one of flavour B plus 
one mode of each flavour antisymmetrised in a singlet, 

,„(6)AB c [A B] c D E F fA a \ 

yjv i ~ £CDEF »1 »t <«f »f Mf »f . l^.oj 



The previous expressions are symbolic and the products of mf's in (|4.4|)-( [4~q ) correspond 
to different combinations of the modes rj A , £ aA , v A and v Au in the various entries of the 
ADHM matrix for each field. The structure of the terms with more fermion modes in the 
multiplet can be determined analogously. The iterative solution of the field equation to 
construct the first few terms in the multiplet was carried out explicitly in [31]. 

From the above equations we can deduce the form of the component fields in the 
decomposition relevant for the BMN limit. For the scalars in (|2.4|) we have 



y(2) [1 4] 7(2) [2 3] 

,,(2)1,3 [1 3] [2 4] ,,,(2)2,4 [1 2] [3 4] Q » 

~ m£ -+- m£ m£ , tp^ ' ~ m £ m£ -f- m£ «i£ . (^•■^J 



For the fermions in ( p. 7] ) and (2.9) we have respectively 

V>~ (1)a ~ m f, ~ i M+m i) a i a = 1,4 o = 2,3 (4.10) 

and 

^ (3) ~e«jcrc«M«? , ^" (3)d ~ (Af-e) & BCJJ»f»f»F , = 1,4,6 = 2,3, 

(4.11) 

so that 

ft 1 ' 1 -^, V~ (1)4 ~«f, ^ +(1)4 ~ m ?, (4.12) 

„/,+(3)2 13 4 „/,+(3)3 12 4 „7;-(3) 2 2 3 4 ,I-(3)3 1 2 3 ( a -1 o\ 
tp y ' ~ m £ m £m £ , 1p y ' ~ m £ m £ m £ , w ~ m £m £ m £ , T/J v ' ~m£m£m£.(4.1oJ 

The terms of higher order are easily deduced from the previous general discussion. Notice 
that we can assign U(l) charge +^ to the fermion modes m\ and m\ and charge —5 to the 
modes mj and mijL 

The dependence on the superconformal modes, r) A and £ aA , can be obtained using 
supersymmetry without solving the field equations. These modes are associated with su- 
perconformal symmetries broken in the instanton background. The corresponding terms 
in the jY=A supermultiplet can thus be generated acting with the broken Poincare and 
special super symmetries, Q a A and SaA, on the classical instanton solution for the gauge 
potential. In the case of SU(2) gauge group there are no additional fermion modes and the 
complete solution can be obtained in this way. In general, however, the dependence on the 
v A and v Au modes can be determined only by solving the equations of motion. 

It is useful to discuss the derivation of the dependence on the superconformal modes 
using supersymmetry since it allows us to clarify how different combinations of fermion 
modes appear in various operators. 

Substituting A^ = A 1 ^ in the supersymmetry transformation of X A gives Aa A > which 
is linear in rj A and £ aA and solves the corresponding field equation. Replacing A^ by X^ A 
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in the variation of ip AB generates the solution if^ AB for the scalar. The iteration of this 
procedure gives rise to A^, then to the correction A^ to the gauge field and so on. 

In the examples studied in [15-17,27] the superconformal modes always appear in the 
expressions of gauge-invariant composite operators in the combination 

= - (x - *o)Xa ^1 • (4.14) 

In general, however, the dependence on the fermion superconformal modes, even in gauge- 
invariant operators, is not only through this combination and instead the moduli r\ A and 
£ aA appear explicitly. This can be understood analysing the form of the Poincare and 
special supersymmetry variations of the fields. Under a combination of the broken super- 



symmetries, r] aA Q a A + w e have 

SA li = (rj A + <T-x^ A )a li X A (4.15) 

5X A = F^(rj A + a-xi A ) + [p AB , <p BC ](r) C + ™f) (4-16) 

5<p AB = \ A ( V B + a-x£ B )-(A^ B) (4.17) 

5\ A = H)(pab(j] B + a-xl B ) + Cp AB l B , (4.18) 



which shows that, whereas the variations of A^, X A and tp involve the combination 
C A , the superconformal variation of X A contains an extra term. Therefore the profiles of 
operators involving X\ in general depend separately on r\ A and £ A . Since further application 
of the broken supersymmetries generates new terms in the solution for the elementary fields, 
it follows that not only operators containing X&A, but also those in which any elementary 
field contain a non-minimal number of fermion modes (e.g. A^\ A« , ip^ AB ) will depend 
on r\ A and £ A not only via ( A . This observation will play an important role in the case 
of two impurity BMN operators. As will be shown in the next section, a naive counting 
of zero-modes including only terms with the minimal number of fermion modes in each 
field would lead to conclude that these operators have vanishing two-point functions in the 
instanton background. We will, however, argue that the inclusion of the term 
the solution is needed in order to compute the leading non-zero instanton contributions to 
these two-point functions. 

5. Two-point functions of BMN operators 

In this section we analyse instanton contributions to two-point functions of the BMN 
operators described in section |2.2| . Using the results of the previous sections we shall 
determine which operators have non-zero two-point functions in the instanton background 
and the dependence of the instanton induced anomalous dimensions on the parameters, 
g YM , N and J as well as the integers corresponding to the mode numbers in the dual string 
states. Zero and one impurity operators are protected, therefore their two-point functions 
are not renormalised and receive no instanton contribution. We shall therefore discuss two 
and four impurity operators. Operators with an odd number of impurities are expected 
not to receive instanton contributions. The analysis of two impurity operators in the next 
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subsection will be rather qualitative because the leading non-zero contribution to their 
two-point functions involves the six-fermion term in the scalar solution which is not known 
explicitly. The four impurity case which is fully under control will be discussed in the 
following subsection. 

5.1 Two impurity operators 

At the two impurity level we focus on the bosonic SO(4)c singlet operator ( |2.17| ) in the 9 
of SO(4)r. Since this sector contains only one operator there is no problem of mixing and 
the anomalous dimension of the operator ^}g. n can be read directly from the coefficient of 

the two-point function ( ^jg n { x i ) ^-j\ m ( x 2 )) ■ 

As usual it is convenient to compute this two-point function for a particular choice 
of components, rather than work in a manifestly SO(4)r covariant way. Therefore we 
consider 

G 9 ( Xl ,x 2 ) = (^ 13 >(xi)^ 13 >(z 2 ))i ns t, (5.1) 

so that there is no trace to subtract. 
The component is 

J 

^{13} = 1 y e ^ P n/j ^ ( Z j-p^ Z p^\ _ Tr (z J "V 4 zV 4 )] (5.2) 

and the conjugate operator is 

J 

V e - 2nipn/J [Tr (Z J ~ V 3 ^V 13 ) " Tr (Z J - V 4 ZV 4 )] • 



(5-3) 

The semi-classical approximation in the one-instanton sector for (|5.1|) gives 



-4N „4N J2mr ,. 4 



G 9 (x 1 ,x 2 ) = ( *l^™ e _ 2) , / d^d 4 x J] dVd^d^Vd^W*- 13 (5.4) 



e 9 ym p &). n J {xi;x Q ,p,r],£,i>,v) &_j.' m {x 2 ;x , p,rj,£,v,v) 



In order to have a non-zero contribution to this two-point function the classical profiles of 
the two operators must contain, when combined, the sixteen fermion modes corresponding 
to the broken supersymmetries. 

It is easy to verify that substituting for each scalar field in (|5.1| ) the leading order 
solution, (f( 2 \ does not allow to soak up all the superconformal modes. Using for each 
scalar field the bilinear term ( f4.9| ) in the solution we find that €?^ 13 ^ contains the following 
combinations of fermion modes 

[1 4]\ J ( [1 3] [1 3] . [2 4] [2 4]\ (v - ^ 



4 The subscripts indicating the SO(4)_r representation and the U(l) charge will be omitted except in 
situations where this may cause confusion. 
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Similarly, the conjugate operator &^ ^ contains 

(5.6) 



.[2 3]\ J fjl 31 [1 3] , 12 4] [2 4] 



Tfl £ 1tl £ P I fft £ fft £ fft £ m £ "1~ fft £ £ £ W £ 



The argument given at the end of section |] shows that in these traces, involving only 
and (p( 2 > AB , the superconformal modes always appear in the combination ( A of ( |4,14 ). This 



can be easily verified using the explicit expression for the scalar ADHM matrices given in 
d A . 8| ) . Then, because of the condition (( A (x)) 3 = satisfied by the Weyl spinors ( A , in 
order to soak up the sixteen superconformal modes in the two-point function ( |5.4j ) each of 
the operators should contain two factors of (, A for each flavour. In other words the sixteen 
superconformal modes should appear in the two-point function in the form 

f[[c A M 2 [c A M 2 . (5.7) 

A=l 



Examining the combinations ( |5.5| ) and (|5.q ) it is clear that this is not possible, 
cannot soak up the required superconformal modes of flavour 2 and 3, since it does not 
contain two factors of both £ 2 and £ 3 , while cannot soak up all the superconformal 

modes of flavour 1 and 4, since it does not contain two factors of both C 1 and £ 4 . This 
simple analysis of the flavour structure of the superconformal modes in the classical profiles 
of the operators shows that the two-point function ( |5.1[ ) vanishes at leading order in the 
instanton background. This argument does not rely on the way the remaining J — 2 (v A v B ) 
bilinears are distributed in the two operators. According to the discussion in section || the 
leading contribution in the large N limit would come from terms in which all the (u A u B ) 
bilinears are antisymmetrised. However, since the above argument is based only on the 
analysis of the superconformal modes the conclusion that the leading g YM contribution to 



the two-point function (5T) vanishes is valid at all orders in 1/N. 

In order to saturate the integrations over the superconformal modes in ( |5.4| ) we need 
to use for some of the scalar fields the solution containing six fermionic modes, ip^ AB . 
Inspecting the combinations ([5l]) and ( |5.6| ) found at leading order and recalling ( |4.8| ) it 
is easy to verify that it is sufficient to consider one (p^ AB (or and Z^ respectively) 
insertion in each operator. These are the leading order contributions, the insertion of more 
six-fermion scalars leads to contributions of higher order in g YM since in this case more 
v A v B bilinears appear. 

Recalling the form of ip^ AB given in ( t4.8|) we find that the combinations of fermionic 
modes in Tr (Z J - p ip AB Z p cp CD ) and Tr (Z J ~ p ip AB Z p ip CD ) are respectively 



A' B' C D' ( [1 4]\ J ( [A B] [C 
EA'B'C D' m f m f m f m f I m f m f I I m f m f m f m 



f ] ) (5.8) 



and 



A' B> C D> ( [2 3]V 7 / [A B] [C D]\ , . 

EA'B'C'D' m { m f m f m f I m f m f I l m f m f m f m f I • l^-^J 

Notice that here the superconformal modes do not necessarily enter via ( |4.14j ), since we are 
using the term with six fermions in the solution for one of the fields in each operator. More 
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precisely the structure of the supersymmetry transformations ( 4.15|) -( f4.18|) shows that both 
traces contain one single £ mode which is not part of a £. This is crucial in order to get 
a non-zero result from the moduli space integration, because it allows, for two flavours, to 
distribute three fermionic superconformal modes in one operator and one in the other. 
The resulting non-zero contribution to the two-point function is 



G 9 (x 1 ,x 2 ) 



4 2ttit 
JiV 3 / 2 



A=l 



+ 



d A x Q dpp 2J 5 /Oi, x 2 ;x ,p) 

(c 1 ) 2 ^ (c 3 ) 2 (c 4 ) 2 £ _1 



d 5 o (n 



23 ^ J -W 



c 1 (c 2 ) 2 (c 3 ) 2 (c 4 ) 2 ? 



(a*) 



(^[(O'CO'COVf]^)}, (5.io) 



where the (v A v B ) bilinears have been rewritten in terms of Q AB, s as described in section 
H The overall powers of g YM and N come from the normalisation of the operators, the 
moduli space integration measure and the {vv)§ bilinears, see equations ( |3.16| ) and ( 3,17| ). 
The fermion superconformal modes are saturated and the corresponding integration is non- 
zero. In ( p,10| ) the dependence on the bosonic moduli has been collected in the function 
f(xi,X2]Xo,p), which can only be computed knowing the explicit form of the solution 
^(6)AB wn j cn we h a ve not determined. The exact form of the solution is also needed in 
order to compute the overall coefficient and, in particular, the dependence on J. More 
details of the derivation of ( 5.10) ) as well as of the evaluation of the moduli space integrals 
are given in appendix [B]. 

The final result for the two-point function is of the form 



G Q (x 1 ,x 2 ) 



j3 e 2niT 



iV 3 / 2 { Xl - x 2 ) 2 ( J + 2 ) 



(5.11) 



where / is a logarithmically divergent integral, to be regulated e.g. by dimensional regu- 
larisation of the xq integral. The logarithmic divergence is due to the bosonic integrations 
over xq and p, as can be verified by dimensional analysis. The presence of this divergence 
signals an instanton contribution to the anomalous dimension of the operator 0g 

As already observed there is only one operator in the representation 9 of SO(4)# and 
thus there is no mixing to resolve and the present analysis directly determines the instanton 
correction to the scaling dimension. We thus find that the instanton induced anomalous 



dimension of 



?{13} 



behaves as 



o 4 J 3 - 

„,mst y YM 



-+iB 



(g 2 f 2 (A') z e-5* 



+i6 



(5.12) 



This is in agreement with the non-perturbative correction to the mass of the dual string 
state computed in [14]. In particular, the anomalous dimension (5.12) is independent of the 
parameter n corresponding to the mode number of the plane-wave string state. Apart from 
the exponential factor characteristic of instanton effects, ( 5.12j ) contains an additional factor 
of (A') 2 . This is due to the inclusion of six- fermion scalars which give rise to additional 
{vv)q bilinears, each of which brings one more power of g YM . As will be shown in the 
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next subsection in the case of four impurity SO(4)# singlets it is sufficient to consider the 
bilinear solution for all the scalars and as a consequence we shall find a leading contribution 
of order ( 52 ) 7/2 e- 8 ^/ 92A '. 

Contributions in which some of the vv bilinears are in the 10 of SU(4) give rise to 
subleading corrections which are suppressed by powers of 1/N. 

Another class of contributions to ( |5.1| ) which are suppressed in the large N limit 
are those in which pairs of scalars are contracted. In these terms the analysis of the 
superconformal modes is unaltered and in order to soak them up it is again necessary to 
use the solution ip(^ AB for two of the scalars. Two of the scalars which were previously 
replaced by vv bilinears are now contracted and do not contain any fermion modes. Hence 
the integration over the moduli space produces one less power of <? 2 M A ■ However, with 
the normalisation we are using the propagator is proportional to g 2 M , so that in conclusion 
the contribution of these terms is down by 1/N with respect to ( |5.11| ) because there is no 
power of N associated with the contraction. 

A careful analysis of both types of 1/N corrections shows that they give a contribution 
to the anomalous dimension of the operator of order (fl^) 9 ^ 2 ^') 2 - These are the leading 
terms in a power series in g2- In general the corrections to the semi-classical approximation 
in the BMN limit can be reorganised into a double series in 52 and A'. 

Operators in different sectors can be studied along the same lines. However, super- 
conformal invariance implies that all the two impurity operators have the same anomalous 
dimension [28] and thus the above result can be extended to two impurity operators in all 
the other sectors with no further calculations required. 

Arguments similar to those discussed here, showing the vanishing of the leading one- 
instanton contribution to the two-point function Gg(xi, x%), have been used to prove various 
non-renormalisation properties in [27,32,33]. In view of the results we found for 0$, one can 
expect that some of the non-renormalisation results of these papers may not be extended 
to higher orders in the coupling. 

5.2 Four impurity operators 

The calculation of two-point functions of four impurity operators is more involved than the 
corresponding calculation in the two impurity case from the point of view of the combi- 
natorial analysis. However, at the four impurity level, in the case of SO(4)_r singlets, the 
leading instanton contributions do not involve the six fermion solution for the scalar fields. 
A non-zero result is obtained using only the bilinear solution, which is known explicitly and 



given in (A.8|), in computing the classical profiles of the operators. Therefore we can anal- 



yse in a quantitative way the semi-classical contributions to the two-point functions. The 
fact that non-zero correlation functions of singlet operators are obtained using the minimal 
number of fermion modes for each field also implies that in this case a contribution to the 
matrix of anomalous dimensions arises at leading order in the instanton background. As we 
shall see these operators have instanton induced anomalous dimension of order (g2) 7 ^ 2 e 2mr . 
Another difference with respect to the two impurity case studied in the previous section is 
that two-point functions of four impurity operators depend explicitly on the integers dual 
to the string mode numbers. We shall discuss in detail a SO(4)/jxSO(4)c singlet with four 
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scalar impurities and show that the behaviour of its two-point functions is in remarkable 
agreement with the corresponding string calculation of [14]. Other singlet operators can 
be analysed in a similar fashion. Operators in other sectors will be shown to receive con- 
tribution only at higher order in A'. This result follows simply from the analysis of fermion 
zero modes and is also in agreement with the string theory prediction. 



5.2.1 e-singlet operator 

In this subsection we present the calculation of the one-instanton contribution to the two- 
point function of one particular SO(4)# singlet. More details are provided in appendix |B]. 
We focus on the four scalar impurity operator in which the SO(4)# indices are contracted 
via an e-tensor, 

J 

/f _ £i 3 kl 2T7i[(n 1 +n2+n 3 )q+(n2+n 3 )r+n 3 s]/J 

'l;J;ni,n2,n 3 — , — - 7^ c 

j3 ( % \ q,r,s=0 



q+r+s<J 



;Tr (z J -^ r+8 ^(p i Z q tfPZ r ^Z 8 (A . (5.13) 



The string state in the plane-wave background which is naturally identified as being dual 
to this operator is of the form 

Sijki al ni al n2 a k : n3 a l __ {ni+n2 _ n:j) \0) h , (5.14) 

where \0)h is the BMN ground state and the contraction runs over values of the indices 
in one of the two SO(4) factors. L>-instanton contributions to the renormalisation of the 
mass of this state were computed in [14]. We shall return to the comparison with the string 
results at the end of this section. Notice, however, that the state (5.14) is antisymmetric 



under the exchange of the two left-moving or right-moving modes. The operator ( 5.13p 
on the other hand has no definite symmetry under permutations of the parameters m, n 2 
and 713. Therefore in order to construct a gauge theory operator that can be identified 



with ( |53g ) it will be necessary to explicitly antisymmetrise (|5,13|). This point will prove 



crucial when comparing instanton corrections to the scaling dimension of 6\ to D-instanton 
induced corrections to the mass of the string state. 
We are interested in the two-point function 



G 1 (xi,X 2 ;ni,n 2 ,n 3 ; mi, m 2 , m S ) = (£i ; , M ,n 2! n3 i x l) ^ljmi.ma.msO^inst 

- [ dpd 4 z ft dVd^d^Vd^VA 13 (5-15) 

/* J A 1 



-AN 4N 2niT 

(N - l)l{N - 2)! 



A=l 



■2 



Xe 169 YM p &l;n u n2,n 3 {xi;X , p, Tj,^, V, u) 01; mu m 2 , m3 (x 2 : X , p, V, V, v) . 

As usual the semi-classical approximation requires the calculation of the classical profiles 
of G\ and G\ in the instanton background. 
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Summing over the SO (4) indices in ( 5.13| ) and using the relations (|2.4j ) we find that 
the operator G\ contains the independent traces 



+Tr (zV 2 ^V 3 ^V 4 ^V 4 ) - Tr {Z^ 12 Z^ u Z r ^ 2A Z s V 13 ) 
+Tr (zV^V^V^V 3 ) + Tr {Z P V 12 Z«<p 34 Z r <p 13 Z s ^ 24 ) 
-Tr (ZV 2 ^V 3 ^V 4 ^V 4 ) " Tr {Z^ 12 Z«<p 24 Z r <p 13 Z s <p 34 ) , (5.16) 

where p = J — (q + r + s), plus three other groups of six traces obtained by cyclic permuta- 
tions of the indices on the impurities in ( [5.16 ), The conjugate operator, 6\, contains the 



same terms, but with the Z's replaced by Z } s. 

It is straightforward to verify that these traces, when evaluated in the instanton back- 
ground, contain the correct combination of fermions required to soak up the superconformal 
modes in a two-point function and that this can be achieved using only the bilinear solution 
for all the scalars. In this case all the r\ A and £ A modes in the gauge invariant traces are 
combined into £ s. In order to give rise to a non-zero two-point function in the one instan- 
ton sector both operators should then contain the combination n^=i (C^) 2 - To achieve 
this in each trace in ( 5.1 Cj ) the four impurities must provide two superconformal modes of 



flavours 2 and 3, whereas the modes of flavour 1 and 4 can be taken from the impurities 
or from the Z J s. Similarly in the case of G\ the superconformal modes of flavour 1 and 4 
come necessarily from the impurities and those of flavour 2 and 3 can be provided by the 
impurities or by the Z's. As in the two impurity case studied in the previous section the 
leading contribution is obtained taking all the remaining modes in (pv)e bilinears. In all 
the traces appearing in both G\ and G\ the impurities contain two fermion modes of each 
flavour. The combination of fermion modes entering into all the terms in G\ is 

(«f) J+2 («f) 2 (»f) 2 (»f) J+2 , (5-17) 
whereas all the terms in the expansion of G\ contain 

( m f) 2 ( m f) J+2 ( m f) J+2 ( m f) 2 • ( 5 T8) 

The leading contribution to the two-point function Gi(xi, X2) in the semi-classical ap- 
proximation arises from terms in the profiles of the operators containing the following 
combinations of fermion modes 

^1 ^ (c 1 ) 2 (c 2 ) 2 (c 3 ) 2 (c 4 ) 2 

^-(c 1 ) a (c a ) a (c 8 ) a (O a (^) J - ( 5 - 19 ) 

As previously observed these combinations can be obtained in many different ways cor- 
responding to the choice of which field, Z or 93, provides each of the £'s of flavour 1 and 
4 in & . An equal number of different contributions arises from the ways of distributing 
the £'s of flavour 2 and 3 among the impurities or the Z's in 6 . In order to simplify 
the discussion of the associated combinatorics it is convenient to introduce the following 
notation. We denote by (p AB a scalar solution in which only the v and v modes are kept 
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and all the superconformal modes are set to zero; scalars containing only bilinears in the 
superconformal modes are indicated by (p~ AB ; the symbol (p AB is used for scalar profiles in 
which only mixed terms, C,v or Qv, are included, 

<p AB = l p AB (x,x o ,p;u,P;r ] = ^=0) (5.20) 
~ab _ ^b (x> p . J] £. v = u = 0) (5 21) 

The same notation is also used for Z ~ 93 14 and Z ~ <p 2S . 

We are only interested in contributions to the two-point function Gi(x\, X2) which 
survive in the BMN limit, N — > 00, J — > 00, with J 2 /N fixed. The leading large N con- 
tributions are those in which the combinations ( |5.19| ) are selected, i.e. the superconformal 
modes are soaked up and all the remaining fields are replaced by (Pv)6 bilinears. Within 
this class of terms the dominant ones in the large J limit are those in which as many 
superconformal modes as possible are extracted from the Z's and Z's, because there is 
roughly a multiplicity factor of J associated with the choice of each Z or Z providing one 
such mode. We first discuss these leading terms and we will then show that these are the 
only non- vanishing contributions in the BMN limit. 

As shown by the previous preliminary analysis, in the operator G\ the modes £ 2 and 
C 3 necessarily come from the impurities and thus the leading large J terms arise from 
traces in which we take the two C} and the two £ 4 modes from four distinct Z's. Similar 
considerations apply to the G\ operator with the role of the flavours (1,4) and (2,3) 
exchanged. Using the notation introduced in ( 5.2C| )- (|5.22 ) this means that we consider 



traces in which all four impurities are (p AB matrices and we choose four Z's to be Z matrices, 
with all the others being Z's. There is a total of 35 different traces of this type for each of 
the 6x4 terms in the operator ( |5.13D and a similar counting applies to its conjugate. The 35 
traces correspond to the inequivalent ways of choosing the four Z's from the four groups of 
Z's in (PTI3D . For the generic term in the operator, Tr (Z p ip AlBl Z q tp A2B2 Z r (p AsB:i Z s ip A4B4 ) , 
with p = J — (q + r + s), we need to consider 

Tr (^Z pl ZZ P2 ZZ P3 ZZ P4 ZZ P5 ^ AlBl ZV A2B2 Z r ^ A3B3 Z s ^ 4B4 ) 
Tr (Z P1 ZZ V2 ZZ Vi ZZ VA (p AlBl Z qi ZZ q2 q(p MB2 Z r (p HBz Z s (p MB ^ 

Tr (z p ^ AlBl ZV A2B2 Z r <^ A3B3 Z Sl ZZ S2 ZZ S3 ZZ S4 ZZ S5 ^ 4B4 ^ , (5.23) 



where in the first trace YliPi = P ~ 4 = J— (q + r + s + A), in the second YliPi = P — 3 
and Yli Qi = 1 ~ 1 an d so on until the last sum where ^ = s — 4. The ellipsis in ( 5.23| ) 



refers to other combinations in which the four Z's are gradually moved to the right. All 
these traces can be evaluated using the ADHM solution for the scalars given in ( |A.8| ) and 
selecting for each factor the matrix elements containing the appropriate fermion bilinears. 
The calculation is rather involved. As explained in appendix [B] it can be carried out most 
efficiently defining a more general trace from which the 35 distinct traces ( |5.23 ) can be 
obtained for different choices of indices. 
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In order to compute the relevant part of the profile of G\ we need to sum the contribu- 
tions of the traces ( 5.23Q corresponding to the 6x4 choices of indices, (Ai, Bi), i = 1, . . . , 4, 
on the impurities. A key feature of all these traces is that they do not depend on the 
way the Z's are grouped, i.e. they do not depend on the exponents, (pi, . . . ,ps, q, r, s), 
(pi, . . . ,P4, q\, q%, r, s) etc. in ( |5.23| ), but only on the ordering of the four Z's with respect 



to the four impurities, (p 



■AiBi 



1, . . . , 4. This is a consequence of the structure of the 



ADHM matrices and the restrictions imposed by the ADHM constraints. Keeping only 
the terms with two £'s of each flavour all the traces in produce expressions which after 
simple Fierz rearrangements can be brought to the form 



[( Xl - x ) 2 + p 2 ]J+* 



(C 1 ) 2 (C 2 ) 2 (C 3 ) 2 (C 4 ) 2 On). 



(5.24) 



Similarly all the contributions from the traces in G\ containing the required eight super- 
conformal modes can be put in the form 



\(x 2 - X ) 2 + p 2 } 



J+8 



(C 1 ) 2 (C 2 ) 2 (C 3 ) 2 (C 4 ) 2 ](X2) ■ 



(5.25) 



Each set of indices (Ai, Bi), i = 1, . . . , 4 on the impurities in each of the 35 traces leads to 
a contribution of the form ( [5.24| )-( |5.25| ) with a different numerical coefficient. 

The fact that the result of all the traces can be reduced to the above expressions 



implies that when substituting into the definition of the operator (5.13) and its conjugate 
a common factor ( |5.24| ) or, respectively, ( |5.25 ) can be taken out of the traces. Associated 
with each of the 35 types of traces there are, however, multiplicity factors which make the 
sums in the definition of the operator non-trivial. For instance in the last trace in ( |5.23j ) 
there are s choices for the first Z among the Z's, (s — 1) choices for the second Z, (s — 2) 
for the third and (s — 3) for the fourth. After substituting into the definition ( [5.13D and 
factoring out the moduli dependence in the form ( |5.24 ), the contribution of the last trace 



in ( 5.23 ) involves the sums 



J 

E • 

q,r,s=0 
q+r+s<J 



2ni[(ni+n2+n3)q+(n2+n3)r+n3s]/ J 



s(s - l)(s - 2)(s - 3) , 



(5.26) 



with a numerical coefficient resulting from the contributions of the 6x4 permutations of 
indices of the impurities. Repeating the same analysis for all the traces means combining 
a huge number of terms which makes the calculation extremely laborious. Completely 
analogous steps go into the calculation of the profile of the conjugate operator. 

After lengthy algebraic manipulations and the use of the formalism described in section 
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|, the semi-classical result for the two-point function ( |5.15| ) takes the form 

e 2nir f^X dp p J+8 p J+8 



G 1 (xi,x 2 ) 



j3 N 7/2 J p 5 [( Xl _ X0 )2 +/O 2]J+8[( X2 _ X0 )2 +/0 2]J+8 
4 4 

A=l B=l 



X 



ridVd 2 ^n[(c B ) 2 H [(c b ) 2 m 

A=l B=l 
Jd 5 n (n l4 ) J (n 23 ) J K(n 1 ,n 2 ,n 3 ;J)K(m 1 ,m 2 ,m 3 ;J), (5.27) 



where following the discussion in section |3| the {yv)§ bilinears have been expressed in terms 
of the angular variables £l AB . In ( [5.27 ) overall numerical coefficients have been omitted. 



The J and iV dependence in the prefactor in ( |5.27| ) is obtained combining the normalisation 
of the operators, the contribution of the measure on the instanton moduli space and the 
factors of g YM y/N associated with {vv)§ bilinears. The origin of the various factors which 
determine the dependence on the parameters g YM , N and J will be summarised shortly. 



The expression (5.27) contains integrations over the bosonic moduli, xq and p, the sixteen 
superconformal fermion modes and the five-sphere coordinates Q AB . The dependence on 
the integers nj, uii, i = 1,2,3, dual to the mode numbers of the corresponding string state 
is contained in the functions K{n\, n 2 , n 3 ; J) and K{m\, m 2 , m 3 ; J). These are given by 



the sum of 35 terms, 

.35 



K(m,n 2 ,n 3 ; J) = ^ c a S a (ni, n 2 , n 3 ; J) , (5.28) 



a=l 



where the symbols S a indicate 35 different sums over the indices q, r, s in which the sum- 
mands are given by the phase factor exp{27ri[(ni +n 2 + n 3 )q + {n 2 + n 3 )r + n 3 s]/J} times 
the multiplicity factors associated with the different distributions of Z's in each case. The 
numerical coefficients c a are obtained combining the contributions of the different permu- 



tations of indices on the impurities for each of the 35 terms. See appendix B.2 for more 
details. 

In the large J limit the leading order contribution to the sums S a can be obtained using 
a continuum approximation by setting x = q/J,y = r/J,z = s/J, so that x,y,z £ [0,1]. 
For instance the sum ( |5,26 ) is approximated as 



J 

e 2ni[(n 1 +n 2 +n 3 )q+(n2+n 3 )r+n 3 s]/ J _ -g^ _ _ 3) 

q,r,s=0 
q+r+s<J 

—> J 7 ^ dx ^ X dy ^ X V d z z ^i(^+^+n 3 )x+(n 2 +n 3 )y+n 4 z] ^ 



Jo Jo 



which shows that it behaves as J for large J. We shall denote by n(ni, n 2 , n 3 ) the function 
of the mode numbers arising from these sums/integrals after extracting a factor of J , 

K(ni,n 2 ,n 3 ; J) = J 7 K,(ni,n 2 ,n 3 ) . (5.30) 
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The two-point function is then 

Jll e 27rir fd A x dp 



G 1 (x 1 ,x 2 ) 



AT7/2 



-,J+S 



.J+8 



p 5 [ (xi _ Xo) 2 +/0 2 ] J + 8[ (a . 2 _ a . o) 2 + p 2 ] 



21 J+8 



iidVd 2 ^n[(o 2 H [(o 2 (- 2 



x J d 5 Sl (£1 14 ) J (0 23 ) J K(ni,n 2 ,n 3 ; J)«(mi, m 2 , m 3 ; J) . (5.31) 

Unlike the case of two impurity operators discussed in the previous subsection, here the 
dependence on the instanton moduli is known explicitly and we can compute the associated 
integrations. More details are given in appendix |B.2j . The integration over the five-sphere 
in ( 5.31 ) is a special case of the general integral ( |3.19| ) and gives 



i s s = Jd 5 n (n u ) J (n 23 ) J 



TT 



(J + l)(J + 2) ' 



(5.32) 



The integration over the superconformal modes is also straightforward. It does not depend 
on iV or J. For each flavour the result is 



dV 2 £ 



(C) 2 (X!) (C) 2 (X 2 ) =-(x!-X 2 ) 2 , 



(5.33) 



so that the fermionic integrals contribute a factor of (x\ — x 2 ) 8 . The integration over the 
bosonic part of the moduli space must be treated carefully since it is logarithmically diver- 
gent as expected in the presence of a contribution to the matrix of anomalous dimensions. 
The integrals need to be regulated for instance by dimensional regularisation of the x$ 
integral and can then be computed using standard techniques, e.g. introducing Feynman 
parameters. The result is 



d xq dp 



.J+8 



.J+8 



p 5 [(xi - X ) 2 + pty+S [{x 2 ~ X ) 2 + p 2 ] 
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(5.34) 



12) 



The 1/e pole is the manifestation of a logarithmic divergence in dimensional regularisation. 
The contribution ( 5.34] ) behaves as 1/J 2 in the large J limit. 

Putting together all the contributions the dependence on the parameters, g YM , N and 
J, in the correlation function can be summarised as follows 
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(5.35) 
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The final result for the two-point function is thus, up to a numerical coefficient, 

2 ue , \ / N 1 

(4.) 



G 1 (x 1 ,x 2 ) = (c/2) 7/2 e ^ ir+te K(ni,n 2 ,n 3 )K(m 1 ,m 2 ,m 3 ) ^ 2 J+4 log (A 2 xj 2 ) . (5.36) 



where the scale A appears as a consequence of the 1/e divergence. It has no observable 
effect. The physical information contained in the two-point function is in the contribution 
to the matrix of anomalous dimensions which is read from the coefficient in ( |5.36| ) and does 
not depend on A. 

The result is expressed in terms of the double scaling parameters A' and g 2 . Note that, 
unlike the two-point functions of two impurity operators ( |5,36j ) is independent of A' apart 
from the dependence in the exponential instanton weight. 



The non-perturbative mass correction computed in [14] for the state (5.14), in terms 
of the same gauge theory parameters, is of the form 

5m ~ ( 52 ) 7/2 e-^^ 9 -^ . (5.37) 

(nin2 r 

so that the A' and g 2 dependence is in agreement with the gauge theory calculation. The 
mode number dependence in ( 5.37| ) is remarkably simple and a special feature of the string 



result, which is a direct consequence of the structure of the .D-instanton boundary state, is 
that it is non-vanishing only if the mode numbers in both the incoming and the outgoing 
states are pairwise equal. The only states which couple to the D-instanton are of the form 

|0) fc . (5.38) 



On the other hand in the gauge theory result ( [5.36 ) obtained for the operator ( 5.13| ) 



the mode number dependence is contained in K(ni,n 2 ,n 3 ) and K(mi,m 2 ,ms), which are 
extremely complicated rational functions of their arguments. In particular the condition 
that the integers rij be equal in pairs does not seem to be required. 

However, as observed after Q5.14 ) in order to correctly match the properties of the dual 



string state, the operator ( 5.13; ) must be explicitly antisymmetrised under the exchange of 



pairs of mode numbers. This antisymmetrisation induces dramatic simplifications. Work- 
ing with the correctly antisymmetrised operators the result for the two-point function is 



Gl{Xl > X2) = ¥^ (n 1 n 2 )(m 1 m 2 ) (^ 2 ) J+4 6 ( ^ ( } 



if the mode numbers in each of the two operators are equal in pairs and vanishes otherwise. 
In ( |5.39 ) we have reinstated all the numerical coefficients coming from the profiles of the 



operators, the combinatorics previous described and the moduli space measure. This result 
is in perfect agreement with the string result (5.37) of [14]. It is worth stressing that the 



simplification found after the antisymmetrisation is extraordinary given the complexity of 
the function «(ni, n 2 , 713). Moreover the condition of pairwise equal mode numbers which 
is also imposed in this way is far from obvious and highly non-trivial from the point of 
view of the gauge theory calculation. 
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As we have seen, in the two-point function computed in the semi-classical approxima- 
tion the mode number dependence factorises. A consequence of this is that in Gx(xi,X2) 
the two independent mode numbers in ff\, n% and do not have to equal those in mi 
and rri2- This appears to contradict energy conservation in the dual string amplitude, which 
requires the mode numbers of the outgoing state to match one to one those of the incoming 
state. However, the fact that the condition m,j = rij, % = 1, 2 does not arise is an effect of 
the semi-classical approximation. This is valid in the A' — > limit which corresponds to 
the m — » oo limit in the plane-wave string theory (where m is the mass parameter entering 
the string action). In this strict limit energy conservation in a two-point string amplitude 
only requires that the number of oscillators in the incoming and outgoing states be equal, 
with no constraint on the associated mode numbers. Therefore ( 5.39| ) is indeed in agree- 



ment with the string theory result. On the other hand the instanton corrections discussed 
here should be considered as subleading corrections on top of the perturbative effects. The 
condition m; = rij on the operators in a two-point function is already imposed at the per- 
turbative level and should therefore be assumed when computing instanton contributions 
in the semi-classical approximation. 

The calculation presented here is not sufficient to determine the actual instanton in- 
duced anomalous dimension of the operator G\. This requires the diagonalisation of the 
matrix of anomalous dimensions of which we have not computed all the entries. Other 
entries are determined by the corresponding two-point functions whose calculation follows 
the same steps described here and results in expressions similar to ( |5.39| ). From this we 
can conclude that the behaviour of the leading instanton contribution to the anomalous 
dimensions of singlet operators is 



_2 



7i nst ~ (92) 7 ' 2 ^^T^ ■ ( 5 - 4 °) 

(nin 2 y 

As a further test of the result we find that the two point function vanishes in the limit of 
zero mode numbers. The function n{ni, 712,^3) is identically zero when n\ = n-i = = 0. 
This is the expected behaviour because in this limit the operator is expected to become 
protected and no corrections to its free theory two-point functions should arise. The string 
theory counterpart of this result is the decoupling of the supergravity modes (dual to the 
protected operators with {m = 0}), which was also verified in [14]. 

In the previous analysis we have considered only a class of contributions in which in 
each operator as many superconformal modes as possible were taken from the Z's. It is 
easy to verify that these are the only relevant terms at leading order in the BMN limit. All 
the other types of traces are suppressed and vanish in the J — > 00 limit. As an example 
consider a contribution to the profile of G\ in which the £ 2 and £ 3 modes as well as one of 
either the C 1 or £ 4 modes are taken from the impurities. Instead of the last trace in ( 5.23f ) 
we would then consider traces of the type 

Tr (^Z p y MBl Z q (p A2B2 Z r (p MBi Z Sl ZZ S2 Z , (5.41) 

where the first impurity contains two superconformal modes and thus only three Z's are 



needed. An analysis similar to that carried out for the traces (5.23) can be repeated in this 
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case and one finds that associated with such a trace there is sum of the form 
J 1 

e ^[( n l+n2+n 3 )q+{n 2 +n 3 )r+n 3 s]/J J_ s ^ g _ -^^ s _ 2^ (5-42) 

q,r,s=0 
q+r+s<J 

which behaves as J 6 in the large J limit. The remaining J, N and g YM dependence in the 
correlation function is unmodified and thus the combined behaviour of this contribution 



can be read from ( 5.35 ) replacing the last factor on the first line by (J 6 ) 2 leading to 

, mT J 5 e 2 ^( 52 ) 7 /2 



c 



N 7/2 J2 



which vanishes in the BMN limit. Similar arguments can be repeated for all the contribu- 
tions other than those leading to ( 5.3SQ , which is therefore the complete leading instanton 



contribution to this singlet two-point function in the BMN limit. We shall briefly comment 
on corrections to this result of higher order in A' and Q2 in the discussion section. 

5.2.2 Other four impurity singlets 

There are two other independent four impurity singlet operators involving four scalar im- 
purities. They correspond to the two inequivalent ways of contracting the SO(4)# indices 
with Kronecker delta's, 

1 J 

j>W0 _ 'y y p 2m[(n\+n2+n 3 )q+(n2+n 3 )r+n 3 s\/ J 



1; J ;ni ,U2 ,n 3 



J 3 



ym \ q,r,s=0 



I 1V1 I 



q+r+s<J 

xTr (z J - ( ^ r+a ^(p i Z q tp i Z r ipjZ'tp') , (5.43) 



1 J 

f( d 2) _ ^ e 2iril(rii+n2+ri 3 )q+(n2+n 3 )r+n 3 s]/ J 



' 73 ( YM \ q,r,s=0 



q+r+s<J 

xTr (z J -(« +r+s yZVZV^V) • (5-44) 



The calculation of instanton contributions to the two-point functions of these operators 
proceeds in complete analogy with the discussion in the previous subsection. These op- 
erators are expected to receive contributions of the same type as the e-sing let ( pU3D and 



to mix with the latter much in the same way as was found in the string theory analysis 
of [14]. 

Other singlet operators can be constructed using all the other combinations of impuri- 
ties in table In all these cases the analysis of fermion zero modes shows that a non-zero 
contribution to the corresponding two-point functions can arise at the same leading order as 



( 5.40 ). Operators containing D^Z insertions correspond to string states involving bosonic 
oscillators which are vectors of the second SO(4). Operators containing the ip~ a and ip~^ a 
fermions are dual to states created by the oscillators. The calculation of two-point 
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functions of all these operators is similar to that described in the previous section with the 
additional technical complication that in the presence of covariant derivatives the solution 
Ay for the gauge potential is needed and for operators containing fermions the solution 
Xa is needed. 

As observed in the case of the operator ( |5.13| ), two-point functions in the semi-classical 
approximation factorise, with the two operators being related only by the five-sphere inte- 
gration. Because of this property mixing is expected among all the operators which receive 
instanton contributions. 

In the ,yV=4 theory it is in principle possible to construct a large number of other 
operators which potentially mix with those considered here, being SO(4)/jxSO(4)c singlets 
with A — J = 4. These involve A — J = 2 impurities and thus do not correspond to new 
independent states having vanishing two-point functions in free theory. However, it is 
known that the inclusion of such operators is needed in perturbation theory to properly 
resolve the mixing beyond the zeroth order approximation in the 52 expansion. Since 
instanton effects are exponentially suppressed in gi one should in principle expect these 
operators to be relevant at leading order in the instanton background. This is, however, 
not the case. The combinatorial analysis involved in computing the classical profiles of the 
operators shows that those containing A — J = 2 impurities are suppressed in the large J 
limit. 



5.2.3 Operators in other sectors 

As observed in section |2.2.2 the spectrum of four impurity BMN operators is rather rich. 



Instanton contributions to the anomalous dimensions of operators in other sectors can be 
studied with the same methods used for the singlets. D-instanton induced amplitudes for 
string states in the plane wave background dual to non-singlet operators are suppressed 
with respect to those in the singlet sector. Hence string theory predicts that the leading 
instanton contributions to the anomalous dimensions of non-singlet four impurity operators 
should be suppressed with respect to (5.40). More precisely the string prediction is that 



the leading non-zero contributions should arise at order e 27TiT (g 2 ) 7/2 {\') 2 - We shall not 
discuss in detail the calculation of two-point functions needed to verify this prediction, but 
we present here an argument indicating that the gauge theory result is indeed in agreement 
with string theory. We focus on an operator with four scalar impurities which is a singlet 
of SO(4) c and belongs to the 3+ 3~ of SO(4) fl , 

1 J 

jfM _ 1 e 2Tri[(n 1 +n 2 +n 3 )q+(n 2 +n 3 )r+n 3 s}/J 

3+®3-;J;ni,n 2 ,n 3 j J+4 

T3 ( YM \ q,r,s=Q 

' ' ^ > q+ r + s<J 

xTr(z J -^ +r+s ^i P k Z q <p k Z r ^ i Z s <p>A . (5.45) 



The study of other non-singlet operators is completely analogous. Considering for con- 
creteness the component i = 1, j = 2 in ( 5.45| ), we find that using for all the scalars the 



bilinear solution the combinations of fermion modes contained in the classical profiles of 
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the operator and its conjugate are respectively 

(<»?) 2 ( 

2\* ( 3\ / 4\ J + 2 
'■{) \ m f) \ m f) 



fj \ m i) \ m i) \ m i) 



/ is J+2 , 2^3 
+ [mf) \m ( 



+ ( m f 



n J+2 



3 \3 / 4n J+2 



+ (« f 



1\ J+l / 2\ 2 I 3\ 2 / 4\^+3 



m f) ( m f) ( m f)' 



and 



+ 



1\3 

1\2 
'fj 



2 x J+2 
'fj 

2 n J+3 
'fj 



3^ J+2 
! fj 

! fj 



4) 

4\2 



4) 

, / IN 2 / 2\ J +f / 3^ J + 3 / 4\ 2 
, / 1\ I 2\J+ 2 ( 3\ J + 2 { 4\3 



(5.46) 
(5.47) 
(5.48) 
(5.49) 

(5.50) 
(5.51) 
(5.52) 
(5.53) 



This shows that the integrations over the superconformal modes in the two-point function 
can be saturated, selecting terms containing ( 5.46j ) or Q5.49 ) in G^ 2 ^ and terms containing 

12 1. However, with these choices the resulting five-sphere integrals 



(I5T5TD and flpl) in 



vanish. For instance combining ( 5.46 ) and ( 5.51 ) leads to the following moduli space 
integrals 
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(5.54) 



The integration over the five-sphere arising from the second line of ( |5.54 ) vanishes because 
the multiplicity of the flavours 1 and 2 exceeds that of the flavours 3 and 4. 

In order to soak up the superconformal modes while avoiding the obstruction from 
the five sphere integral it is necessary to include a six-fermion term in each operator. In 
this way the combinations of modes in the two operators are the same as in ( |5.46| )-( |5.53l ) 
with the addition of one mode of each flavour. The same arguments given in section ^4] in 
connection with two impurity operators can be repeated here and for instance combining 
( 5.46 ) and ( 5.53j ) we get moduli space integrations of the type 



iVe [(c x ) 2 c 2 (c 3 ) 2 (c 4 ) 2 e _1 ] (xi ) [c 1 (c 2 ) 2 (c 3 ) 2 (c 4 ) 2 ?] 

x J d^W 4 ^ 2 ^ (VVl) J+1 (vK^) J+1 (VVl) (VM) . (5.55) 

Just as in the two impurity case the resulting non-vanishing contribution to the two-point 
function is suppressed by a factor of (A') 2 due to the additional {i>v)§ bilinears in ( p. 55 ). In 
conclusion the analysis of fermion zero modes confirms that the leading non-zero instanton 
contribution to the anomalous dimensions of four impurity operators in the 3 + and 3~ 
representations behaves as 
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(5.56) 
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in agreement with the string prediction of [14]. 

Other sectors can be analysed in a similar fashion and we find that all non-singlet 
operators receive leading non-zero contributions of the same order as ( 5.56 ). 



6. Discussion and conclusions 

This paper has considered one-instanton contributions to two-point correlation functions 
of gauge invariant operators in the BMN sector of the jV=A supersymmetric Yang-Mills 
theory. These determine the leading instanton contributions to the anomalous dimensions 
of operators dual to physical string states in the maximally supersymmetric plane-wave 
background obtained as Penrose limit of AdSs x S 5 . The basic message is that we find 
striking agreement between these instanton effects in the gauge theory and those of the 
plane- wave string theory calculated in [14]. 

We focused on operators with two and four scalar impurities. The four impurity case, 
although more involved, is fully under control, whereas the two impurity case presents 
subtleties due to the fact the leading semi-classical approximation vanishes and the first 
non-zero contribution arises at higher order. We have explicitly computed a two-point func- 
tion of four impurity operators which are SO(4)^jxSO(4)c* singlets. Our analysis shows that 
instanton induced contributions to the anomalous dimensions of operators in this sector 
behave as l/(nin2) 2 exp (—8ir 2 /g2\' + id) g^ 2 , where A' and <?2 are the effective coupling 
constant and genus counting parameter in the BMN limit and n\ and rii correspond to 
the mode numbers of the dual string state. The result is in perfect agreement with the 
D-instanton correction to the mass matrix elements of the corresponding states in the 
plane- wave string theory which was computed in [14]. Even without directly matching the 
numerical values of the anomalous dimensions and the string mass renormalisation, the 
agreement with the string calculation appears highly non-trivial. The correct dependence 
on the parameters A' and g2 is obtained by combining contributions arising from the in- 
tegrations over the instanton moduli space and various combinatorial factors. Even more 
impressively, the mode number dependence found in [14] is reproduced after spectacular 
cancellations. 

In the case of two impurity operators the leading instanton correction vanishes. The 
first non-zero contribution is awkward to calculate completely, but with mild assump- 
tions we showed that it has the form exp (—8n 2 /g2\' + iQ) g\ A' 2 and does not depend 
on the single mode number characterising the dual string state. This behaviour is also 
in agreement with the results of [14], although the subtleties presented by the gauge the- 
ory calculation did not arise on the string side. Four impurity operators in sectors other 
than the singlet have the same A' and g% dependence as two impurity operators, again in 
agreement with the string prediction of [14]. 

Our results provide a significant new test of the duality proposed in [1]. The fact 
that non-perturbative contributions obey BMN scaling, i.e. can be re-expressed in terms 
of the effective parameters A' and g% , strongly supports the conjecture that this property 
should hold at all orders. This can be further tested by analysing subleading effects in 
the instanton background. A class of higher order contributions can easily be obtained 
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from the calculations presented in this paper, relaxing the requirement that all the fermion 
modes of type v and v be combined in bilinears. As already observed, replacing 

a bilinear with a (vv)\q leads to a suppression by a factor of l/y/N, see ( 3.17] ). 



The profile of BMN operators of the type that we considered contains J (Pz^)'s after the 
superconformal modes have been soaked up. Hence each factor of 1/yN coming from the 
replacement of a (Pu)q by a (i'v)io is associated with a factor of J corresponding to the 

1/2 

number of choices for the (z/^)io bilinear, resulting in a suppression by g 2 • Moreover in 
order to get a non-zero result from the five sphere integration an even number of (z^)io 
is required. Therefore contributions of this type with an increasing number of (z^)io 
insertions give rise to subleading corrections which form a series in integer powers of g2- 
More complicated terms with the same behaviour correspond to contributions in which 
pairs of fields are contracted between the two-operators. We can also identify a class of 
subleading corrections suppressed by powers of A'. These are generated by including in 
the profile of the operators a number of fermion modes greater than the minimal number 
required by the moduli space integration. For instance in the case of the four impurity 



singlets that we studied in section 5.2.1 this is achieved by including one six- fermion scalar 



in one of the two operators. In this case the calculation is analogous to that of section 



5.1| for the two impurity case and we can argue that the resulting contribution should be 
suppressed by a factor of A'. Including more six- fermion terms gives rise to higher powers 
of A'. Although these arguments are rather qualitative they indicate how the perturbative 
series of corrections to the semi-classical one-instanton contributions can be reorganised 
into a double series in A' and g2- 

Some rather striking general properties of the D-instanton induced corrections to the 
string mass spectrum of [14] can be immediately deduced from the structure of the D- 
instanton boundary state in the plane-wave background, whereas the corresponding effects 
in the gauge theory are far from obvious. In fact, the string theory results suggest a 
number of extensions and generalisations of the gauge theory results. For example, one 
generic feature of the string calculation is that only states with an even number of non-zero 
mode insertions receive D-instanton corrections. Zero mode oscillators can appear in odd 
numbers with the condition that they be contracted into a SO(4)xSO(4) scalar between 
the incoming and outgoing states. The simplest example in which these properties can 
be verified involves five impurity operators and the calculation of the necessary two-point 
function in the gauge theory is extremely complicated 5 . In general, contributions to 
operators with a larger (even) number of impurities are expected to be non-zero at leading 
order in the instanton background. However, the complexity of the combinatorics involved 
in such calculations grows rapidly with the number of impurities. 

Another peculiarity observed in the string theory calculation [14] is that the D- 
instanton contribution to the masses of certain states with a large number of fermionic 
non-zero mode excitations involves large powers of the mass parameter m. When ex- 
pressed in terms of gauge theory parameters this corresponds to large inverse powers of 
A'. As observed in [14] the behaviour of these mass corrections is not pathological in the 

5 Three impurity operators present technical difficulties similar to those encountered in the two impurity 
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A' — ► limit, because the inverse powers of A' are accompanied by the instanton factor 
exp ( y —87r 2 /g2\')- From the point of view of the gauge theory this result is particularly in- 
triguing not only because of the unusual coupling constant dependence that the anomalous 
dimensions of the dual operators should display, but also because there are no other known 
examples of operators in ^V=A SYM whose anomalous dimension receives instanton but 
not perturbative corrections. We will study this particular class of BMN operators in a 
future publication [25]. 

In the original formulation of the AdS/CFT duality, relating jV=A SYM to type IIB 
string theory in AdSs x S 5 , the effects of multi-instantons in the large- A^ limit of the 
gauge theory and of multi D-instantons in string theory were shown to be in remarkable 
agreement [16]. Clearly it would be of interest to generalise the present work from the one- 
instanton sector to the multi-instanton sector. However, such a generalisation is technically 
very challenging both on the string and on the gauge side. 

Instanton effects have been studied in a number of different supersymmetric field theo- 
ries [34-40] in the context of the AdS/CFT correspondence, and agreement has been found 
between string and gauge theory. The example of the JY=2 Sp(A r ) superconformal field 
theory studied in [34,37,38] is particularly interesting in connection with our work because 
in this case the analogue of the BMN limit has been studied in [41]. In this case the duality 
involves a theory of open and closed strings in a plane-wave background and the dual gauge 
theory has a rich spectrum of gauge-invariant operators and possesses a Higgs branch [42]. 
It would be interesting to study instanton effects in the BMN sector of this theory. 

In a conformal field theory, the problem of computing the scaling dimensions of gauge 
invariant operators can be reformulated as an eigenvalue problem for the dilation operator 
of the theory. At the perturbative level this observation leads to a very efficient approach 
to the calculation of anomalous dimensions in jV=A SYM [43]. Some comments about the 
possibility of extending this approach to non-perturbative sectors were made in [27], but 
there has been no further progress in this direction. A remarkable consequence of recasting 
the problem of computing anomalous dimensions as an eigenvalue problem for the dilation 
operator is the emergence of connections with integrable systems. In the planar limit the 
dilation operator can be related to the hamiltonian of an integrable spin chain, leading 
to the possibility of applying techniques such as the Bethe ansatz to the computation of 
anomalous dimensions [44], see [45] for a review and references. The integrability structure 
observed in the ^V=A Yang-Mills theory appears, however, to be spoiled by the inclusion of 
non-planar contributions. Therefore instanton effects, which are exponentially suppressed 
in the large- iV limit, are unlikely to be relevant in connection with integrability. 

Instantons play a special role in the jV=A theory in connection with the SL(2,Z) S- 
duality symmetry, which transforms the complex coupling so that r — ► (with a, b,c,d£ 
% satisfying ad — be = 1) thereby mixing perturbative and non-perturbative effects. In 
the conformal phase invariance of the theory requires that the full spectrum of scaling 
dimensions be invariant. The anomalous dimensions are therefore naturally expressed as 
functions 7(7-, r). Similarly D-instantons are instrumental in the implementation of S- 
duality in type IIB string theory. Their role is well understood at the level of the effective 
action for the supergravity states, but little is known at the level of the massive string 
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excitations. As in the SYM case, invariance of the theory requires that the complete 
spectrum be invariant. In general SL(2,Z) transformations relate operators of small and 
large dimension, just as in string theory they relate fundamental strings to D-strings, which 
have large masses of order l/g s , in the limit of weak string coupling, g s <C 1. It would be 
interesting to understand how 5-duality is realised in type IIB string theory in the plane- 
wave background. A corresponding symmetry should exist in the BMN sector of jV=A 
SYM and the instanton effects which we studied in this paper should be important in its 
implementation. 
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A. Useful formulae 

This appendix contains some definitions and formulae used in the paper. The Clebsch- 
Gordan coefficients T> AB (£>f B ) projecting the product of two 4's (4's) onto the 6 are 
defined as 

^ AB = ^ AB ,^ A + B 3 ) = (r ] a AB , i f] AB ) 

^t B = ^ab^a + b) = (-rf B ,itf B ), (A.l) 
where a = 1, 2, 3 and the 't Hooft symbols jf AB and fj AB are 

Vab = Vab = £ aAB , A, B = 1,2,3, 
ri a A 4 = fj a 4A = S a A, 

Vab = -V%a , Vab = -fj a BA ■ (A.2) 

In some situations the JV=\ formulation proves very useful. The JV=\ decomposition of 
the t /f / =4 supermultiplet consists of three chiral multiplets and one vector multiplet and 
under this decomposition only a SU(3)xU(l) subgroup of the SU(4) R-symmetry group is 
manifest. The six scalars are combined into three complex fields, (j) 1 , I = 1, 2, 3, according 
to 

4=i=(^-^+ 3 ) . (A.3) 

The complex scalars (p 1 and (f)\ transform respectively in the 3i and 3_i of SU(3)xU(l) 
(where the subscript indicates the U(l) charge). The fermions in the chiral multiplets are 

Vd = A^, # = A?, 7 = 1,2,3, (A.4) 
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transforming in the 33/2 and 3_3/2- The fourth fermion and the vector form the JV=1 
vector multiplet, {X a = ,A^}, and are SU(3)xU(l) singlets. 

Using these definitions we find the following relations among the scalars in the different 
formulations 



-2 V2, 



^ = ^ (V 4 + ^ 23 ) , ^ 
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(V 3 + v? 24 ), 



2 

i\/2 



(^ 12 +^ 34 ) 



(A.5) 



and 



^ = 2^ , <p 2 = 2^ , 



J.3 _ n 34 



2ip , 



4 = 2<^ 23 , 4 



-2v9 13 , 4 = 2 <P 12 



(A.6) 



In the ADHM formalism the expressions for the t /( / =4 elementary fields in the back- 
ground of an instanton are conveniently given as [N + 2] x [N + 2] matrices. In particular, 
the two-fermion solution for the scalar field (p AB in the one instanton sector can be written 
in the block-form 



(l)AB\ v 



(2)AB\ 7 
/ u; 



where u, v = 1,2, . . . , N , a, (3 = 1,2 and 



(1)AB\ v 



^{y 2 [-m &B if-e% 



(4)AB^ 7 J ' 



(A.7) 



4(y 2 + p 2 

+{y 2 + p 2 ) [-4(£M - + - ^ Bv )} 



+ y 
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'/3 



(2)AB\ 7 



A; 4(y 2 + p 2 ) 2 



16y Vd^V - f*V B ) + 4(y 2 + p 2 )^ 



W u ,f 



ie y &s (vF^ A - ntrf B ) + \f^v & \v Au v B - } 



(3)A_bA « 



^2 {p 2 [l6y^(f B § - f^f)™^ - 4ypa(l 



/3; 4(y 2 + p 2 ) 

- 16fo?& - r^f )™ d; * + 4(t/|^ - ^)] } 



-Ad 



A-Bv\ 



+ l y —±/ L q } {v Ar v B -v Br 4) 



2 p 2 

B. Instanton induced two-point functions of BMN operators 



(Ai 



In this appendix we present some details of the calculations of one-instanton contributions 
to the two-point functions of BMN operators discussed in section 13. 
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B.l Two-impurity operator in the 9 of SO(4)^ 



As shown in section the leading semi-classical contribution to the two-point functions 
of two impurity operators in the 9 of SO(4)_r vanishes because the superconformal modes 
cannot be soaked up. A non-zero result is obtained including for one scalar field in each 
operator the six fermion solution. 



In the case of the component considered in section 5.1 the terms in the two-point 
function which contain the correct combination of fermion modes to give a non vanishing 
contribution are 

_ , > 1 / 2iripn\ ( 2iriqm\ 

G 9 (x 1 ,x 2 )= ^ cos ( j~ ) cos ( J ) 

jfWlj P , q =o V J ' \ J J 

(Tr [(Z^V^V 3 ) (xi)] Tr [(Z J ~VW 4 ) (x 2 )]) 

(Tr [(Z J "V^V 4 ) (xi)] Tr [(Z J ^ 13 Z^ 13 ) (x 2 )])l . (B.l) 



The other terms vanish in the instanton background either because they do not contain all 
the required superconformal modes or because of the integration over the five-sphere. For 
instance if one considers (Tr [[Z J ~ p (p 13 Z p ip 13 ) (xi)] Tr [(Z J ~ q (p 13 Z q (p 13 ) (x 2 )]) it is easy 
to verify that the superconformal modes can be soaked up, but the resulting five-sphere 
integral vanishes because among the remaining fermion modes different flavours appear 
with different multiplicities. 

Let us consider the terms in flB.ip where in each trace one scalar is understood to be 
replaced with the six-fermion solution. In the first expectation value the two traces contain 
respectively the following combinations of fermion modes 

[m { ) {m t ) [tn { ) (m f J 

( m f) 1 ( m f) ( m f) ( m f) 3 • ( B - 2 ) 

Using the fact that each trace contains one £ mode not part of a £ we can soak up the 
superconformal modes selecting the following combinations of fermion modes in the two 
traces 

^V 2 (fi\ 2 (r^ 2 ^ /V,[i,/]V _1 (rM 



J-l 



Tr(z j -v 3 ^v 3 ) - (c 1 ) c 2 (c 3 ) (c 4 ) 

Tr (Z^V^V 4 ) -> C 1 (C 2 ) 2 (C 3 ) 2 (C 4 ) 2 f (v [2 v 3] Y~ L (v [2 v 4] ) • (B.3) 



Similarly the two traces in the second term in (B.l) contain respectively 

(m}) J+1 (m 2 ) 3 (m^) 1 (mj) J+3 (B.4) 

and 

( m f) 3 ( m f) J+1 ( m f) J+S i" 1 ^ 1 ( B - 5 ) 
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and we need to consider 



Tr (Z J ~V 4 ZV 4 ) - (C 1 ) 2 (C 2 ) 2 C 3 (C 4 ) 2 e 4 ^ '" Jl " 4 
Tr(Z J -V 3 ^V 3 ) -> (C 1 )* (0* (C 3 )Yf (^ 3] 



J-i 



J-i 



(B.6) 



These expressions contain the correct combinations of fermion superconformal modes such 
that the corresponding integration is non-zero. The two terms in (|B.1|) give rise to 



II dVd 2 f{[(cf( 2 (C 3 ) 2 (t 4 ) 2 e](xi) [C 1 (C 2 ) 2 (C 3 ) 2 (C 4 ) 2 f ](x 2 ) 

4=1 

+ [(C 1 ) 2 (C 2 ) 2 C 3 (C 4 ) 2 e 4 ] (xi) [(c 1 ) 2 (c 2 ) 2 (c 3 ) 2 C 4 f ] (x 2 )} 

~ (x 1 - x ) ■ (x 2 - x Q ) (xi - x 2 ) 4 . (B.7) 

After re-expressing the (i> A v B ) bilinears in terms of Q AB, s as described in section [3] both 
(B.3) and (|B.6|) lead to the same five-sphere integral, 



d 5 o (o 1 



(0 



This can be evaluated rewriting it as 



'55 



Jr. 



d 6 o (s^) J_1 (sf ^y -1 (4 3 o fc ) (e 24 o 



(B.8) 



(B.9) 



where the symbols are defined in ( AT ). Defining O = £* 4 O l = (O 1 + iO 4 ), SI = 

£ 23 0* = (O 1 - iO 4 ), O = £f O' = (SI 2 + iO 5 ) and O = £ 24 0* = (O 2 - iO 5 ) the integral 
reduces to 



i s b = / d 6 o 5 o 2 - 1 (oo) 



J-i 



00 



dO dO dO dO d 2 7 ft J + 00 + 00-1) (00) J 1 (^00 
where O 7 = (0 3 ,0 6 ). Introducing polar coordinates for the SV directions 

J s6 = 2vr /" dr r dO dO dO dO 5(r 2 + 00 + 00-1) (00) J_1 (00 

dOdOdOdo (oo) J_1 (oo) . 



(B.10) 



7T 



The remaining integrals are straightforward 



(B.ll) 



I S 5 = 2vr 

^2 



dOdO (00) 



j-i 



dzz" 



o 



3 fl 

— I dO dO (1 - 00) 2 (00) J_1 = — / dz z (1 - z 2 ) 
2 Jnn<i ' " 

3 



2,* Z 2J 







7T 



J(J + l)(J + 2) 



(B.12) 
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As observed in section 5A the exact dependence on the bosonic moduli in the two-point 
function Gg(xi, X2) cannot be determined without knowing the six-fermion solution. Di- 
mensional analysis indicates that the final bosonic integrations over position and size of 
the instanton are logarithmically divergent as expected in the presence of an instanton 
contribution to the anomalous dimension of the operator G^ 1 ^ . 

B.2 e-singlet four impurity operator 

In order to compute the profiles of the operator G\ in ( |5.13D and its conjugate, which are 
needed in the calculation of the two-point function ( |5.15| ), we have to evaluate the traces 
( 5.23 ). In the instanton background these are rewritten as traces of [N + 2] -dimensional 
ADHM matrices. To compute these traces more efficiently it is convenient to define the 
[N + 2] x [N + 2] matrix 



v,/3 



(B.13) 



where the notation used is that introduced in ( |5.20| )-( |5.22j ). This has the standard block- 



form of ADHM matrices and the range of the indices here is the same as in ( |A.8[ ) for the 
elementary scalar fields. 

In terms of the matrix U k "j^ % ' j ' D * ((, v, P) all the 35 traces we are interested in can be 
written as 



Tr 



U, 



Ci,Di;C2,£)2, 



Gs,Ds;Ce,De . 



ki,k 2 



(B.14) 

for appropriate choices of the indices Cj, Di and the exponents hi, i = 1, . . . , 8. For example 
the three traces written explicitly in ( |5.23j ) become 

Tr MifAC, u, "Wl^AC, v, y)U^ A ^ u, u)U^ A ^(C, v, P)] , 



Tr Kf/iC, », y)Ultit' Bl (C, m 



(C,v,")u r A ! 



B 3 ;A 4 ,B 4 



and 



A 1 ,B 1 ;A 2 ,B 2 , 



(c, v, v)u a h b ^\q, v, ^'%' 4 (c, * m 



rl,4;A 4 ,B 4 
s 4 ,s 5 



The generic trace ( |B.14| ) is thus the only one that needs to be evaluated. It can be computed 
using the building blocks (A. 8) and the result is 



Tr 



ttCi,Di;C 2 ,D 2 
U ki,k 2 

1 



08 



C 3 ,D 3 ;C 4 ,D 4 



2 3J-8 [( x _ Xo )2 +p 2]J+8 

+ ( ^ Cli U l] 



J -8 



C5,D5,Ce,De , 



r^Di^D 2 \ r^D 3 ^D 4 \ ^D 5 ^D 6 \ (( D 7£ D s)] 



+ [(C D2 C Di ) (C D K D5 ) (C De C D7 ) (C Ds C Dl )] 



y[C 3 ^] 



„[1 V C4] 



+ I 1; ■<'•■ ir 



y[C2^] 



y[Ce^} 



ylCSy*] 



yllyC?] 



+ p 



v \\c 3 \ 



yl^Cr] 



yl*yC 6 ] 



D lCr v 4] 



yllyCs] 



yllyCs} 



+ permutations > , 



(B.15) 
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where the permutations not indicated explicitly correspond to antisymmetrisation in all 
the (Cj, Di) pairs. 

The key feature of ( B.15 ) is that it depends on the set of indices (Ci,Di), but not 



on the exponents fcj. The traces ( 5.23|) require four pairs of (Ci,Di) indices to be (1,4) 
while the remaining four pairs correspond to the (A/., B^) indices carried by the impurities. 
The fact that (B.15) does not depend on the ki's means that the traces ( 5.23j ) do not 
depend on the exponents on the Z's, but only on the relative positions of the Z's with 
respect to the impurities. Therefore when substituting into the definition ( |5.13| ) of the 
operator the traces ( p.!5| ) can be taken out of the sums over the indices q, r, s. After 
substituting the values of the indices corresponding to the various terms in the expansion 
( 5.16 ) and some simple Fierz rearrangements the profile of the operator G\ takes the 



form of a common factor containing the dependence on the bosonic and fermionic moduli, 
multiplying the combination K{n\,n%,Uz\ J) of 35 sums which contain the dependence on 
the mode numbers n\, and 713, see ( 5.28|) . To illustrate more concretely how this works 
let us describe explicitly one particular term. We consider the first trace in ( |5.16| ) and 
compute the contribution of the last type in (|5.23p for this trace. We have to evaluate 



E 



_,27ri[(rti +n 2 +n 3 )q+ (n 2 +n 3 )r+ri3s] / J 



n 2 7V\ J + 4 

%M \ q,r,si,...,ss=0 



"8^~ ) 



q+r+S!-\ \-s 5 <J-4 

siH \-s 5 =s— 4 



X 



E 



z p,-Ki[{n\+n2+n3)q+(n2+n3)r+ri3s\/ J 



a 2 iy\ J+4 
• YM \ q,r,s 1 ,...,s 5 =0 
8tt 2 J 



q+r+s 1 -\ hs 5 <^~4 

siH \-s s =s—4: 

^{ul^U^Ul^Ul^. (B.16) 



Using ( |B.15| ) with the particular choice of indices in the trace in ( |B.16 ) we get (up to a 
numerical constant) 



rjY /rrl,2;l,3 r r2,4;l,4 r rl,4;l,4rrl,4;3,4\ 
V p,q r,si W S2,S3 w S4,ss J 



2 3J+8 [( a . 1 _ XQ )2 + p 2y+8 



2 2 3J + 8 [(xi - x ) 2 + p 2 ' 



J+8 



{^ ] Y '[(C 1 ) 2 (C 2 ) 2 (C 3 ) 2 (C 4 ) 2 ] , (B.17) 



where the last line has been obtained using simple Fierz rearrangements on the £'s. As 
anticipated the trace is independent of the exponents, q, r, Sj. Equation ( B.16| ) then 
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becomes 

3 1 1 



x ^ e 2^[(ni+n 2 +n3) 9 +(n 2 +n 3 )r-+r l 3 S ]/ji. s ( s _ 1 )( s _2)( s _3) ) (B.18) 
5,r,s=0 

where we have used the fact that there is no dependence on the single exponents, si,. . . ,s§, 
so that the result is independent of the way the four Z are distributed among the last s 
Z's. This leads to the factor ^s(s — l)(s — 2)(s — 3) which is a multiplicity coefficient 
associated with the number of ways of picking four identical Z's out of s Z's. Equation 
QB.18 ) illustrates the factorisation of the result into two terms, the first line containing the 
dependence on the instanton moduli and the second line containing the dependence on the 
mode numbers. 

The function K(nx,n2,n 3 ; J) takes the form 

35 

K(n 1 ,n 2 ,n 3 ; J) = ^ c a S a (ni, n 2 , n 3 ; J) , (B.19) 

0=1 

where each of the S a (ni, n 2 , n 3 ; J) is a sum similar to the second line of (B.18) with dif- 
ferent summand corresponding to the different multiplicity factors associated with the 
distributions of Z's in the traces ( 5.23j ). Table || summarises the contributions to ( p. IS ). 

Using the coefficients given in table [||, and noting that the phase factor factorises, 
( B.19| ) can be written as 

3 3 

K(nx,n 2 ,n 3 ) = -- ^ e 27ri[(n 1+ n2+n 3 ) q +(n 2 +n 3 )r+n 3S ]/J ^ p _ q + r _ s )(p + q + r + s )3 
p,q,r,s=0 
q+r+s+p=J 

oj3 J 

_ ou e 2Tri[(n 1 +n2+n 3 )p+(n 2 +n 3 )q+n 3 j-]/ J _|_ 2r — J) . (B.20) 

p,q,r=0 
p+q+r<J 

The sums in (B.20) can be approximated with integrals in the J — > oo limit, which can 
then be evaluated differentiating a generating function. The relevant generating function 
is given by 

0(CH, 03,03,04) = f dx f X dy \ X V dz^Haix+^y+asz+a^l-x-y-z)] 

Jo Jo Jo 

= -ly . (B.21) 

Thus the above sum requires evaluating 

Si f d d \ 
/(01, 02,03,04) = — ( 2 ^- + 2 q^ ~ 27Ti ) g{ai,a 2 ,a 3 ,a 4 ) . (B.22) 
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Summand 


Coefficient (c a ) 


Summand 


Coefficient (c a ) 


^p(p-l)(p-2)(p-3) 


-18 


m prs(s - 1) 





±p(p-l)(p-2)q 


-9 


±ps(s- l)(s-2) 


9 


±,p(p-l)(p-2)r 


-18 


^g(g-l)(g-2)(g-3) 


18 


±p(p-l)(p-2)s 


-9 


±q(q-l)(q-2)r 


9 


j2T^p{p- l)q(q-l) 





Jjg(g-l)(g-2)« 


18 


l / -i\ 
2\P(P- W 


-9 


jWpQKQ- l)r[r- 1) 





l / 1\ 

2\P{p- i)qs 





2!<7(<7 - 


9 


jzrppip - l)r(r - 1) 


-18 


j2\^q{q- l)s(s- l) 


18 


1 / 1\ 

2\P(P- l)rs 


-9 


igr(r-l)(r-2) 


-9 







2<qr{r - 1)8 





±pq(q-l)(q-2) 


9 


1 / 1\ 

j,qrs(s - l) 


9 


2]Pq(q - i)r 





±qs(s - l){s - 2) 


18 


2<pq{q - i)s 


9 


ir(r-l)(r-2)(r-3) 


-18 


i i i \ 

2;pgr(r - 1) 


—9 


^r(r- l)(r-2)s 


-9 


pqrs 





^r(r - - 1) 





^pqs(s - 1) 


9 


£r S ( S -l)( S -2) 


9 


^pr(r-l)(r-2) 


-18 


^s(s - l)(s - 2)(s - 3) 


18 


ipr(r - l)s 


-9 







Table 3: Contributions to the function K(n\, ri2, n^; J). S a are sums over the indices g,r,s6 [0, J] 
(with the constraint q+r+s < J) in which the summands are those indicated in the table multiplied 
by the phase factor exp (2iri[(ni + ri2 + na)q + (112 + n^)r + n^s]/ J). Here p — J —(q + r + s). The 
c a 's are the combined coefficients taking into account all the terms in the operator. 



As discussed in section 5.2.1 in order to get the correct mode number dependence we need 
to antisymmetrise ( B.22 ) with respect to the exchange of pairs of mode numbers. Therefore 
we need to compute 



"3- 



lim [/(m + n 2 + n 3 , n 2 + n 3 , ra 3 , 0) - f(m + n 2 + n 3 , n\ + n 3 , n 3 , 0) 



-ni 



-/(-ra 3 , -n\ - n 3 , -n\ - n 2 - n 3 , 0) + /(-n 3 , -n 2 - n 3 , -ni - n 2 - n 3 , 0)] 

8ir 2 (nin 2 ) ' ^ ^ 

where only in the case where we impose pairwise equality do we get a non-zero result. In 
conclusion the mode number dependence in the profile of the operator G\ is 

3 J 7 



K(ni,n 2 ,ns; J) 



(B.24) 



87T 2 (nin 2 ) ' 

where the factor of J 7 is the combination of the J 3 in (|B.20 ) and a J 4 arising from the 
conversion of the sums into integrals in the continuum limit. 
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The two-point function Gi(xi,X2) thus becomes 

Jll e 2nir 1 fd*X dp p J+8 p J+8 



G 1 (xi,x 2 



N 1 / 2 (nin 2 )(mim 2 ) J p 5 \(x\ - x ) 2 + p 2 ] J + 8 [(x 2 - x ) 2 + p 2 ] J + 8 

x / fl dVd 2 ^ II [(( B ) 2 (*i)] [(C B ) 2 (a*)] [ d 5 ft (ft 14 ) J (ft 23 ) J . (B.25) 

J A=l B=l J 

The integrations in the second line of ( B.25| ) are straightforward. The integrals over the 
fermion superconformal modes give (x\ — x 2 ) s , see ( |5.33 ). The five sphere integral is similar 
to that encountered in the two impurity case and can be calculated in a similar fashion. 
Proceeding as in ( p.8| )-( p37T2j ) we get 

i s5 = J d 5 n (n l4 ) J (n 23 ) J = J dndnd 4 n r s^n 1 + (nn) J , (b.26) 

where ft = (ft 1 + ift 4 ), ft = (ft 1 - ift 4 ) and ft 7 = (ft 2 , ft 3 , ft 5 , ft 6 ), so that introducing 
spherical coordinates 

I s b = 2vr 2 J dr r dft dft 5(r 2 + ftft - 1) (ftft) J 

= 2tt 2 / dftdft((l-ftft)(ftft) J = — — ^— -. (B.27) 



The integration over the bosonic part of the moduli space must be treated carefully since 
it is logarithmically divergent as expected in the presence of a contribution to the matrix 
of anomalous dimensions. The integrals need to be regulated for instance by dimensional 
regularisation of the xo integral. Introducing Feynman parameters we get 

d 4 x dp p J+8 p J+s 

~p^~ [On - x ) 2 + p 2 } J + 8 [(x 2 - x ) 2 + p 2 } J + 8 



T(2J + 16) 



/ daida 2 5(1 — a\ — a^a^ a 2 + ' 
Jo 

2J+u 



^ X ° P [(xo - a\X\ - a 2 x 2 ) 2 + p 2 + aia 2 x\ 2 ] 2J+1Q 



After dimensional regularisation, 



7h _ _ r(J + 6)r ( j + 8 + e ) 2 _ e i f 1 

[IV + 8)] 2 (xhY^Jo 



a(l - a)] 1+e 



_ i r(j + 6)r(j + 8 + g ) 2 _ e i 
6 [r(J + 8)] 2 (x 2 2 )^ + ,- ^ 

The 1/e pole corresponds to a logarithmic divergence in dimensional regularisation. 
Substituting into (B.25) we finally get 



2 



Cta.*,) ~ fe)" 2 e-i^ + ' (niBa)(mima) ( ^ )J+4 log(A^ 2 ) , (B.30) 
where the exact numerical coefficient was given in ( |5.39| ). 
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